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1. OVERVIEW 


1.1 STUDY 

World wide, continuous measurement of lightning location, 
intensity (or energy) and time during both day and night is to be 
provided by the Lightning Mapper (LITMAP) instrument. This study 
was primarily a technology assessment to determine if the LITMAP 
requirements can be met using existing sensor and electronic 
technologies. The baseline concept discussfd in thi final 
report is a compromise among a number of opposing requirements 
ground resolution versus array size; large field of view 
versus narrow bandpass filter). Th<’ concept does provide 
coverage for more than 80* of the lightning events as based on 
recent above cloud NASA/U2 and satellite lightning measurements. 

1.2 DRIVERS 

The elements making up the Lightning Mapper (LITMAP) sensor 
system are illustrated in Figure 1.1. Those elements enclosed 
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FIGURE 1.1 LIGHTNING MAPPER TECHNICAL ELEMENTS 




within solid lines fall within the bounds of the study statement 
of work (SOW). The primary development risk areas are indicated 
with an asterisk and include focal plane preprocessing, focal 
plane readout and residual data processing. The remaining 
elements enclosed in the dashed lines, (source characteristics, 
optics, telemetry and utilization) did not have to be defined in 
detail as an official part of this study. 

The prime scientific drivers are rapid frame time coupled 
with high spatial resolution and oper tion in a night and high 
background day environment while maintaining a low energy 
detection threshold. These scientific requirements translate 
into the engineering requirements of very high data rate and 
pixel by pixel background subtraction. Additional technology 
drivers identified during the system trades were: focal plane 
array type and size, narrowband filter and a wide field of view. 

1.3 CONCEPT SELECTION 

As with any application, there are multiple approaches that 
can be taken which must then be evaluated against schedule, cost 
and risk factors. A minimum cost, low risk approach for LITMAP 
was selected as the baseline design from which parametric systems 
trades were made. An approach using a single mosaic array could 
be implemented but it would have up to several tens of kilometers 
undetected strips where the mosaic elements are butted together. 
Operationally, slight field of view directional adjustments could 
provide the desired full coverage. An alternate approach for a 
single optical train is to use fioeroptic field splitting. Both 
of these approaches face a large field of view (FOV) versus 
narrow band filter technology problem. This technology problem 
can be solved by either reducing systems performance requirements 
(25A bandwidth rather than 5A) or higher development risk 
(blrefringent instesd of Interference filter). Thus, there is 
latitude to implement alternate variations of the baseline 
concept ss driven by final user considerstions , trades and 
desires. 
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1.*t BASELINE DEFINITION 


The low risk baseline sensor concept discussed in the 
following sections is built around existing imaging CCD arrays 
which require minimal design and development. Since the optics 
aperture is only around ten centimeters, separate telescopes are 
used to minimize the FOV of each array thus permitting the use of 
state-of-the-art narrow band interference filters. The focal 
plane must be partitioned to handle the data rate, but RCA has 
suggested to TRW an approach which leaves the focal plane sensors 
intact and only modifies the output mux, hence minimizing 
development cost and risk. 

1.5 FINAL REPORT CONTENT 

Section 2 summarizes three sets of data: a. LITMAP design 
goals as extracted from the statement of work (SOW); b. baseline 
design parameters and; c. an optimistic parameter set which 
should be attainable with only a slight increase in risk. The 
optimistic parameters set could be pushed even further if a 
development effort is implem nted to combine the best attributes 
of the various technologies and manufacturer's techniques. 
Details defining these parameters are provided in Appendix A. 

Section 3 highlights the system parametric trade 
conclusions. The ' supportive details for this section are 
provided in Appendix A. 

Section A provides the concept selection and baseline 
concept description. The objective was to provide sufficient 
concept definition to permit design trades, t'chnology 
assessments and flight development costing. 

The technology assessment is given in Section 5.0 where 
device availability, cost, schedule and complexity are evaluated. 
Primary emphasis was placed on the focal plane array with 
parallel evaluations for the optics, narrow band filter and 
advanced processing techniques being performed under separate 
IR&D activities 

The flight hardware development and manufacturing cost are 
covered in a separate volume. 
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ORIGINAL PAGE FS 

2. SYSTEMS PARAMETERS OF POOR QUALITY 


2.1 REQUIREMENTS 

Exhibit A entitled "scope of work" (SOW) for the lightning 
mapper (LITMAP) sensor design study provided specification goals 
and lightning characteristics as a general guide for the 
technology assessement of large array sensors and background 
handling techniques. The primary LITMAP program objective is to 
provide measurements of the energy or power of individual 
lightning strokes along with thtir time and location of 
occurence. These measurements are to cover the full disc with a 
90% probability of detection and a false alarm rate of less than 
10%. A summary of these goals and characteristics is provided 
in Table 2.1. 

From preliminary analysis, the primary system drivers are in 
three areas: very high data rate resulting from full time, high 
resolution coverage coupled with a few milliseconds of framing 
time; background clutter during daytime operation due to solar 


TABLE 

2.1 Statement 

of Work 

Specification 


AREA 

ITEM 

VALUE 

URIT 

CATEGORY 

COffflEWTS 

Source Six* 

Souro* Slze-Dlameter 

5 

km 

Characteristic 

typical Cell Size 


Resolution 

5 

km 

Goal 

fully Illualnated Cell 


ritld of View 



Goal 

Full Disc 8 Conus 

Source Power 

Dynaalc Rang* 

10 7 -10 12 

Watts 

Goal 

Total Optical Power 


Power Distribution 

10 7 

Watts 

Characteristic 

90S Greater Than 


Threshold 

10 7 

Watts 

Goal 

SOS Prob. Detection 

Statistic* 

Detection Efficiency 

90 

Percent 

Goal 

of AU Events 



<10 

Percent 

Goal 

of Total Events 


Pulse Duration 

<1 

Nilsee 

Characteristic 


Other 

Wavelength 

8683 

7778 

Angstroms 

Angstroms 

Characteristic 

Characteristic 

MI (1) 
01 U) 


Date Dissemination 

<5 

Minutes 

Goal 

Near Real Time 


Risk-Cost /Schedule 



Goal 


Output 

Limp 

Location 

Intensity 

Tie* 


Goal 

Nap or L, t, T. 

scatter in 

g from clouds 

; and the dynamic range of 

events. These 

factors are interrelated and must 

be weighed 

carefully when 


PRECEDING PAGE BLANK NOT FILMED 
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OF POOR QUALITY 


optimizing the system. The overall objective of the system 
trades and concept definition is to show that daytime lightning 
detection can be accomplished. 

2.2 DESIGN PARAMETERS 

The primary focal point for a systems analysis of an 
optical/detector system performance capability is the signal to 
noise ratio (SNR). This SNR relationship relates the various 
parameters of the detection system through which system 
optimization can be pursued via parametric analysis and 
performance trades. The SNR is defined as follows. 

The number of optical signal generated carriers, N St in one 
sensor element (pixel) is given by: 


N s = * *s 

4 (f/»o) 2 


K Q T? h v 5 • 

E PH0T 


The terms are defined in Table 2.2. Two factors have been 
added to accommodate the statistical lightning event splitting: 
spatial ( v ) and temporal ( $ ). These factors will be Hscussed 
later . 


The noise generated carriers come from several sources 

including those due to the signal (N s ) and the background (Nb)> 
(e.g. solar illiminated clouds). The detector also provides 
several contributions including amplifier (N arop »rms), and dark 
noise (N dn> rms). These sources are assumed to be uncorrelated, 
hence will add in quadrature to give the total rms noise N n of: 


N 

n 


V N s + n b * + 

V^T - V i + c 


V N B V 1 + < N S + N AHP * N dr' /N B 


where the latter form is useful during parametric analysis 
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whenever <<1. The 
background irradiance 


number of carriers 
is given by: 


created 


by 


the 


solar 


"O'b T , 

4 (f/HoV 


AX Alb • 
E PH0T 


TABLE 2.2 Definition of Terms 


OPTIMISTIC 


gnggL 

BUSLINE 

(CHANGES) 

UNITS 

DESCRIPTOR 

Alb 

0.7 


— 

Albedo (scene reflectance) 

A* 

121 x 10 6 


-* 

Lightning Source Aren 

A| 

131 x 10 fl 


■* 

Background Area 

C 

3.00 x 10* 


m/m 

Vacuua Speed of Light 

% 

4 x 10" fl 



Source Radiance (goal) 


3.20 x 10- 19 


J 

Energy per Pbotoo ■ he /A 

r/# 

1 



Effective F/# 

b 

6.63 * 10' 34 


Js 

Planke Conetant 


300 


f/e 2 -sr-ue 

Background (eolar) Irradiance 

i 

0.4 

0.5 


Optical Tranaalaeloa 

Q 

0.35 

0.5 


Detector Quant ue Efficiency 

T| 

4 x 1Q~ 3 


•«e 

Integration Tlee 

X 

M 

5.4 x 10“ 4 


•ec 

Event Width FWHM 

b 

3 x 10~ 5 


■ 

Source laage Six# 

■ 

4 P 

n 

3 x 10~ 5 
0.15 

» 

o 

K 

■ 

Detector Pixel Six# 

Lightning Spectral Efficiency at A 

A 

5.553 x 10“ 7 


■ 

Center Wavelength 

AA 

5 x 10 -4 


um 

Spectral Bandpass 

V 

0.53 


— 

Geometric Split Factor 

C 

(1 - X /3Ti ) 


— 

Temporal Split Factor 

II 

ISO 


e'e 

Readout Noise 

“P 

R 

3.55 x 10 7 


■ 

Sensor Altitude 


The signal to noise ratio is defined as the peak signal to 

the rms noise and is given as: 

SNR = N $ /N n 

V ^ Ng (1 + e) 
which upon rearranging becomes 


- (%’) 



z(r/Mo> 


VT 



1 
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The factors in the SNR expression can be grouped in a number of 
ways, one of which is illustrated above with the segregation 
being based on the physical origin of each term (c.f. Figure 
2.1). Factors 1,2 and part of 6 are source/object related. Item 
3 is primarily optics/filter related; and 4 and 5 are dictated by 
.ien^ur characteristics. Items 6 along with the total pixel count 


FIGURE 2.1 PRIMARY SYSTEM ELEMENTS AND ASSOCIATED PARAMETERS 


S2SrC« 

A. 


"phot 


Ag/pix«l 


rov 

r/# 
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AX 


Array 81m 
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and sensor array partitioning quantifies the post detection 
electronic processing bandwidth. 

Each term in the SNR relationship is discussed and 
quantified in Appendix A and summarized in Table 2.2. With the 
practical constraints defined, performance trades can be made 
with the objective of optimizing the design based on a given set 
of crite-ia (c.f Table 2.1). 

''ince a desired SNR is a priori known, i.e, it is derivable 
fr^ > the desired average false alarm rate (FAR) and event pulse 
cjration (T L ), the chosen approach was to solve the SNR equation 
for the equivalent attainable threshold energy. The parametric 
trades then focused on determining the practical combination of 

l 
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parameters which will satisfy as many of the requirements in 
Table 2.1 as possible. 

The following section summarizes the results of these 
parametric trades. 

2.3 SUMMARY 

Table A.4 summarizes the results of each parameter evaluated 
in Appendix A. This set of data provided the baseline 
configuration input to Table 2.2 as well as the "optimistic” data 
set. Two (often neglected) terms, geometric and time splitting, 
have been included in this analysis. As shown in Section 3 and 
Appendices A and B, these factors can influence the threshold 
performance capability by a factor of two or three, a very 
significant consideration when the total system is pushing the 
technology limits as LITMAP is. 


1 PM WOT rnwgp 


9 



3. SYSTEMS ANALYSIS 
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3.1 DAYTIME OPERATIONS 

As stated in Section 2, the SNR is determined by the FAR 
specification. Thus the approach used for the parametric study 
was that of determining the attainable spectral threshold energy 
for each set of parameters. The SNR relationship can be inverted 
to give the threshold energy, NE s ^ f for a discrete spectral line 
providing : 


he 


st 


SNR 




2(f/No) _ 1 _ 

vr 




PHOT 


I B Alb 


W VTT7 . 

( 1 -^/ 3 ^) 


Section 2 and Appendix A discuss and define the limits for 
each of the terms in this relationship. The above form of the 
expression is most useful when the sensor system noise is 
dominated by solar irradiance for which case e-*0. As amplifier 
noise becomes more prominent, the threshold value is modulated by 
the € term. A practical worse case condition would be when the 
non-solar noise sources equal the solar contribution for which 
€=i and sr 1 + e - 1.41,1. e., a forty-one percent rise in the lower 
threshold limit relative to the condition where solar is the only 
major "noise" contributor. 

Two groups of variables appear in the E s relationship; the 
linear more influential parameters (SNR, f/No.) and the those 
whose influence is somewhat suppressed via the square root 
dependancy. Background/source area (modified by the geometric 
splitting factor) and the time (with time splitting) lie 
somewhere between these two groups. 

3.2 NIGHT OPERATION 

By definition of SNR 
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Day : 


SNRj 




N AMP + N S0L + N s 


Night: SNR 


N 


sn 


/i- 


■W * "s 


The relative day-to-night signal electrons (for constant SNR) is 
given by: 


Vs, 


V 


n amp * N so) 


+ "s / 



sol 


r l + c 


But near threshold N g « N 2 amp + Ngol in day 

and N s << N 2 arnp at night, hence 


Vsn 8 V 


N "aMP + N sol / N AMP 


/ 


If N 


amp 


;VvT/ 


3 then 


Vs, 8 


Y" S0 1 (1nt ^ ,, AMP (nw) 


N 2 amp = N so i ( a probable worse case design condition) 


Nsd/Nsn 3 - 1.41. 


As indicated previously, the nighttime threshold would be about 
41% lower than the daytime. This small improvement is due to the 
large magnitude of the amplifier noise contribution. 

If N 2 amp ^ Ngoi^ 10 then 

n SD/Nsn =>/NsoiVN ara p = 3 

that is, the nighttime threshold would be a factor of three lower 
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than the daytime value. For a value of N go ^ s 500K, the 
amplifier noise would have to be less than 236 electrons rms. 
This noise level appears to be attainable. An upper limit of 500 
e rms has been projected (Hughes, 5 MHz operation, custom device) 
and would give a N sd /Nsn * 1-73, i.e., a 73% improvement for 
nighttime performance over daytime operation. 

If a significant improvement in nighttime operation is 
desired, the non-solar background contributions to the overall 
noise must be minimized and an automatic threshold adjustment 
incorporated into the electronics. 

3.3 ADJUSTABLE THRESHOLD 

The results of the system analysis of available technology 
against LITMAP requirements indicate that the solar background 
will be a severe limitation during the daytime. Reaching the 
goal of 4 u J/m^-sr at the 90% probability of detection with less 
than 10% FAR is very demanding. An array designed with large 
pixel size (hence large aperture) and large full well (FW) (hence 
longer integration time) will drive the system capability in the 
right direction. 

There is another technique which would help optimize the 
threshold setting versus data output rate. The average sensor 
array threshold can be adjusted based on event throughput rate. 
The throughput rate would be set at a value which would be 
compatible with system constraints such as the down-link rate, 
onboard storage capability, etc. As the background energy 
decreases the throughput event rate will decrease for a fixed 
threshold setting. The threshold can be lowered until the event 
rate is back up again (a lower limit may also be desirable). 
When the storm activity increases, the threshold would be 
automatically adjusted to avoid overloading the down-link, etc. 
In this way, maximum sensitivity can be maintained which, under 
some conditions, may not meet the flase alarm (FAR) or 
probability (Pj) specifications but would provide the maximum 
amount of data consistant with the systems constraints. The data 
could be further reduced on the ground using sophisticated 
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analysis techniques. For scientific purposes this may prove to 
be very beneficial in that aspects of the FOV nonuniformities 
could be further suppressed, hence permitting extracting 
information which may have been lost if a fixed threshold had 
been utilized. This type of threshold adjustment is being used 
on a pixel by pixel basis for the newer satellite surveillance 
systems. A combination of fixed minimum plus an adjustable 
threshold level would contribute toward optimization of the 
LITMAP sensor system. 

3.4 PARAMETRIC TRADES 

The primary functional dependency of the detector threshold 
setting on most of the parameters can be deduced from the E s 
relationship given in Section 3.1. Once an absolute baseline 
value is established, effects of other parameters such as F/No., 
A\ , and SNR can be determined by direct scaling (linear, 


TABLE 3.1 Parametric Analysis Summary 
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square root etc.). For parameters such as time and spatial 
splitting, the effects are not as easily visualized, hence 
parametric computations, graphically presented, become very 
useful. Table 3.1 summarizes some of the parametric analyses 
performed and Appendix B contains further discussions. 

3.5 KEY SENSOR PARAMETERS 

Although a number of parameters can be varied (e.g., F/No., 
SNR), the two primary adjustable parameters for the sensor array 
are pixel size and full well capacity. As indicated in Appendix 
A, these two parameters are somewhat related in that the larger 
pixel sizes have larger full well capacities. The relationship 
between these parameters and the minimum resolvable lightning 
threshold setting is given by: 


FW(min) = Bp2(min) s (DR N s (min) ♦ N s 2 (min )/SNR2) 

>- "max signal ma7 noise 

6p(min) = B“1 (DR N s (min) ♦ N s (®in)2/SNR 2 ) 

N s r SNR 2 1/2 ♦ SNR >/sNR2 z 2/4 ♦ N2 amp 

Z = (1 +AAIb t Alb/nE a ) 

B = 1.1 10l5( max ), 0.37 lO^min) 

Figure 3*1 provides *)E S versus pixel size and full well 
requirements for several dynamic ranges. The shaded area for a 
100:1 dynamic range accounts for the variability of existing 
devices. This variation is primarily related to device geometry 
and can be varied by design for a given configuration. The 
dashed lines illustrate that for a given threshold, larger pixel 
sizes (and FW capabilities) are required for larger dynamic 
ranges. Between 30 to 50 micron pixel dimensions with a full 
well (FW) capacity of 10 6 electrons are required to permit a 
source of 4 10“ 6 J/m 2 -sr to be detected from geosynchronous orbit 
using the baseline system specifications in Section 2. Smaller 
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FW capacities result in high thresholds. This figure illustrates 
that the NASA/Galileo 15 micron pixel configuration is not 
compatible with the baseline LITMAP full disc requirements. 
Existing Westinghouse, Texas Instruments (TI) and RCA devices 
come closer to the desired specifications but 3till require 
partitioning to accommodate the short frame time requirement. 



FIGURE 3.1 EASELINE DEVICE PIXEL AND FULL WELL REQUIREMENTS 


3.6 SOLUTIONS TO SYSTEM DRIVERS 

The parametric trades indicate that the system drivers can 
be eased by implementing some or all of the following (full disc, 
geosynchronous orbit assumed): 
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Background Clutter 


o Narrow band optical 
filter 

o Large aperture 

o Larger full well 
capacity 

o Smaller ground 
resolution 

o Background removal 


less background noise: 

gain in SNR : 

■* permits larger 
aperture : 

-*■ larger aperture, 
lower background: 

■+■ eases dynamic range: 


<5A 

Pixel > 20 um 

>5 105 e's 

< 7 km 
frame 

subtraction 


Dynamic Range 

o Larger pixel larger full well: >>20 um 


Data Rate 

o Device selection 

o Partitioning 

o Larger integration 
time 

o Taylored output 
amplifier 


-*■ high clock frequency: 
lower rate per port: 

-*• longer frame time: 

■+ lower readout noise: 


>5 MHz 
>16 ports 
>H msec 

<300 e/pixel 


These factors are discussed in this report 


The degree to which the full potential of the imaging CCD 
technology can be realized depends on a number of factors. The 
main sensor factors are pixel size coupled with total pixel 
count. Other subsystem constraints include wide FOV couoled with 
narrowband filter and large pixel count coupled with short 
framing time. 


As discussed in this report, there are a number of acceptable 
approaches which can be taken to implement the LITMAP Instrument 
with minor modifications to existing devices. Table 3*2 
summarizes some of these configurations. The first entry is the 
baselined approach discussed in this report from which all 
parametric trades were made. It provides the full disc coverage 
with the best performance/cost ratio with the large-FOV/narrow- 
bandwidth-f liter combination being a key driver. 
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TABLE 3.2 Summry of Sensor Configurations 
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Custom imaging CCD arrays with single optical trains are 
covered in the second group of configurations. Although this 
group is within the present technology, the costs and risks are 
higher due to the requirement for larger CCD arrays which are not 
now in production and also due to the requirement for a SOLC 
oiref r ingent filter to handle the large FOV and narrow bandwidth 
requirements. The sensor and filter can cost a factor of 2 to 3 
above the corresponding baseline items. This cost/risk penalty 
may be deemed warrented by the ultimate users. 

The third group contains two approaches which are not deemed 
suitable or practical for the full disc LITMAP application. The 
last group does not provide continuous full coverage but car. 
provide r.ear continuous coverage of specific storm areas using 
unmodified existing devices with a considerable decrease in 
electronics . 

Section B.11 evaluates a few existing commercial imaging CCD 
devices which contain the necessary on focal plane frame storage. 
The results for a 4 millisecond integration time are: 


Device 

uJ/m.2-sr 

Device 

uJ/m.2-sr 

RCA 

9 (x2) 

RCA 

11 (x3) 

Tl 201 

12 (x3) 

Westinghouse 

12 (x3) 

Westinghouse 

14 (x3+) 

TI 201 

17 (x4) 



TI 202 i 
Hughes * 

20 (x5) 


The (x2) refers to a threshold which is about a factor of two 
above the NASA measured 4 uJ/n>2-sr. A (x2) to (x3) would be 
tolerable if the signal dynamic range is 100:1 or greater. 

The large pixel, RCA device gave the best daytime 
sensitivity. Array partitioning and tayloring of the output 
amplifier for high speed operation are needed. The TI devices 
operate at TV rates (> 5 MHz) but array partitioning is still 
required to get the desired framing time. The Westinghouse unit 
is small (100 x 100 x 2) and would require a mosaic of 64 to 100 
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units to handle full disc coverage from geosynchronous orbit. 
This unit would be suitable for a Shuttle or U2 operation. 

The RCA or TI sensor devices would provide the minimum risk 
and cost (< 1M) for implementing the LITMAP instrument using 
multiple optics. 

Special larger imaging array devices can be produced which 
match the LITMAP requirements more closely with an increase in 
development risk and cost (> 2M + risk). For a single optical 
train, the large FOV may require a birefringent narrowband filter 
which may approach a million dollar development cost. 

TI has had experience in producing the larger imaging CCD 
arrays followed by RCA. Both of these companies are preparing 
arrays for the commercial TV market. Hughes, Fairchild and 
West inghouse also produce imaging CCDs. Hughes has provided a 
preliminary estimate to TRW for the development of sensor arrays 
which are compatible with the baseline configuration discussed in 
this report (c.f. fourth entry in Table 3.2). 

3.7 SUMMARY 

The primary conclusion based on a) 90% of the above cloud 
data is above 4 u J/m^-sr and b) the intensity dynamic range 
between 10%-90% is at least 100:1 is that significant lightning 
data can be collected from geosynchronous orbit. Existing 
commercial devices can approach the LITMAP performance 
requirements of sensitivity (pixel size/aperture) and dynamic 
range (full well capacity). The integration and cycle time can 
only be attained within the present technology via securing 
partitioned sensor arrays with the total elements per port on the 
order of 10,000 to 40,000 pixels. This partitioning is quite do- 
able but requires custom devices. 
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4. CONCEPT DESCRIPTION 


4.1 GENERAL 

This section deals with baseline hardware selection and 
concept definition. Trade discussions are briefly repeated to 
support the selection of a given hardware implementation 
approach. The details leading to the selection of the baseline 
set and the parametric trade results are summarized in Sections 2 
and 3 with the bulk of the data given in Appendices A and B 
respectively . 

The previous discussion of lightning characteristics, 
scientific requirements and signal/noise considerations have 
illustrated the challenge surrounding the design of an instrument 
to map lightning from synchronous orbit over the full disk of the 
Earth. The following sections discuss a proposed concept which 
satisfies the lightning mapper objectives. 

The four major considerations are: 

o The optics, including the telescope, bandpass filter, any 
possible image splitting, and the detector at the focal 
plane. In the following section on optics, the issue of 
interference versus birefringent filters is discussed and 
three schemes for partitioning the image among multiple 
detectors will be weighed. 

o The detector. This section is primarily devoted to 
possible CCD formats with alternatives discussed briefly. 
The emphasis will be on creating a device which satisfies 
the requirements with minimal modification to existing 
designs. An argument is made for using multiple sensors 
to cover the full earth disk. 

o Signal processing includes that part of the instrument 
which takes the raw signal from the detector, decides 
when a lightning event has occured and returns a digital 
signal proportional to lightning intensity. The 
discussion will cover the background suppression 
technique, choice of analog or digital background memory, 
and a method of adaptively setting the lightning 
detection threshold to allow improved nighttime 
performance. 

o Data formatting. For every lightning event, a data word 
is created into which is packed the lightning intensity, 
location, and time of occurence. These data are saved in 
a small local memory for telemetry at a low data rate. 
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4.2 OPTICAL CONCEPT 


The LITMAP performance criteria impose signif ics’.t design 
requirements on the optical subsystem. Thus the optimisation of 
the total sensor system requires an interactive trade between the 
optical and sensor constraints. The optical design effort was 
not a part of this NASA study but was evaluated under a separate 
IR&D activity. The following provides an overview of this 
effort . 

4.2.1 Requirements 

The desired field-of-view is a square 17.4 degrees on a side 
covering the entire earth disk. For a subsatellite ground 
resolution of about 13.6 km, an instantaneous angular FOV of 3.8 
x 10” 1 * rad is required. This roughly corresponds to 800 x 800 
elements. 

Use of existing arrays limit the size of the focal plane to 
2.4 cm with a pixel size of 30 um. The effective focal length of 
the system must then be 7.8 cm. Concurrently, the radiometric 
calculations for the system indicate that an entrance pupil of 
15.2 cm is necessary. The resultant system has an F/No of 0.5. 
Although computer simulations indicate that such a design is 
possible for monochromatic light, in practice it would be 
prohibitively expensive. Section 2 defines some of the optics 
parameters which provide a workable compromise between scientific 
requirements and technology constraints. The filter transmission 
bandwidth of the system must be less than 5A to obtain sufficient 
background supression. 

The table below summarizes the above ideal requirements and 
also the baseline specifications presented in Section 2. The 
following sections provide a description of some of the 
alternatives considered for the filter and optics. 
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72 \ 


® 


Baseline 

9.7 cm 

9.7 cm 

0.256 

(14.9 deg) 

3.21x10“^ rad 

5 A 

8683 

4.2.2 Telescope Concepts 

A catadioptic front end telescope (c.f. Figure 4.1) is 
typically used for this type of application. Clear apertures up 
to 15 cm were evaluated via computer. The image from the 
telescope of Figure 4.1 was minified 1:4 using a lens cluster. 
All goals can be met except for the spectral bandwidth. The best 
that an interference filter could do in such systems is 30A 
bandpass which is a result of the large FOV requirements. The 
filter constraints are discussed later. 

An all refractive approach to a wide field telescope was 
also designed (c.f. Figure 4.2) for which 80% of the energy falls 
within a blur diameter corresponding to approximately one half of 
a resolution element. This approach also provides most of the 
desired features except for the filters. Again the F/No of the 
internal telescope combined with the FOV limits the bandpass of 
interference filters to more than 20A. 

To accommodate the filter bandwidth constraint, various 
schemes for splitting the FOV after the primary images were tried 
and are discussed later. In all the splitting schemes the relays 
were mounted on an axis of symmetry for the partial FOV resulting 
in angular shifts of the filter bandpass on the order of 5 A. 


Parameter 


Ideal 


Aperture 

Effective Focal Length 
Total FOV 

Instantaneous FOV 
Spectral Bandwidth 
Spectral Line Center 


15.2 cm 

7.8 cm 

0.303 rad 
(17.4 deg) 

3.78X10- 2 * rad 

< 5 A 

868 ? 
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FIGURE 4.1 CASSEGRAIN 



F/3 OBJECTIVE ♦ 1/4.5 TRANSFER LENS 


FIGURE 4.2 REFRACTIVE 

The simplest approach to the FOV/bandwidth problem is to 
design a refracting telescope with 1.28 rad (7.4 deg) total FOV. 
Then a cluster of four of these 10 cm diameter telescopes would 
cover the entire FOV (c.f. Figure 4.3). Several advantages are 
apparent. First, existing 400 x 400 CCD arrays can be easily 
used at the four independent focal planes. Second, filters may 
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be placed over the apertures or at intermediate surfaces in the 
primary lens cluster and have maximum wavelength shifts of only 5 
A for a required bandpass of about 3 A if the center wavelength 
and filter shape are placed properly. The design of the filter 
will be discussed in the section on filters. This latter 
approach (4-cluster) was selected as the baseline concept. 

4.2.3 Optical Partitioning 

The lightning mapper is a system in which the optical field 
splitting clearly improves performance and the feasibility of 
fabrication.. 

The three approaches selected for inclusion in the study 

were : 

1. Four telescopes, 

2. Knife edge image slicing, and 

3. Fiber optics. 

Each was shown to be completely feasible but possessed differing 
degrees of complexity and risk. 

o The four telescope approach was chosen as the baseline 
design as mentioned previously. The choice was based on 
minimizing the ri3k associated with system development. 



FIGURE 4.3 FOUR CCD'S WITH INDIVIDUAL OPTICS AND SPECTRAL FILTERS. 
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The alternate approaches are as follows. 

o Galileo Electrooptics indicated that a quad ri furcated 
image plane 10 x 10 cm could be fabricated with 20 micron 
spatial resolution. The primary image plane Figure 4.4 
would be divided into four implemented focal planes each 
4.8 cm on a side. Depending on signal processing needs 
each focal plane could be imaged onto an 800 x 800 array. 

A further division into 16 separate channels would allow 
use of 16 600 x 400 arrays to fill the fiber optics area. 
This would achieve twice the ground resolution compared 
to the quad-telescope and hence enhance the threshold 
capability of the system. The alternative to using more 
detector arrays is to reimage the output from the fibers 
onto the 400 x 400 arrays with a 1:4 minification (Figure 
4.5) as in the basic telescope. Only a single collector 
is required with either 4 filters or a single wide F0V 
filter. 

o The image slicer is a knife edge mirror in the focal 
plane of the primary optics (c.f. Figure 4.6). Since 
this only divides the images into two parts additional 
1:1 relay lenses and knife edges are required for the 
final division into 4 fields (Figure 4.7). Each field 
must be minified and imaged as before and utilizing 
either 4 filters or a single wide F0V filter. 



FIGURE 4.4 QUADRI -FURCATED FIBER OPTICS 
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FIGURE 4.6 IMAGE SLICER, TWO WAY 
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FIGURE 4.7 FOUR-MAY IMAGE SLICER CONCEPT 
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4.2.4 Filter 
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The threshold performance of the LIPMAP, when the scene is 
sunlit, pivots on background radiation rejection. Consequently, 
the spectral bandpass of the system must match the line shape of 
the Nitrogen line at 8683 A as closely as possible. 

Unfortunately, with interference filters and birefringent 
filters an increase in the cone angle and/or field-of-view for a 
system increases the effective bandpass of the filter. The ideal 
situation would be to have all the light incident on the filter 
at a single angle. The limitations of the two filter types were 
investigated both theoretically and with potential manufacturers. 

M ulti-c avit y Interferenc e Filt ers The two cavity 
interference filters provide both a narrow spectral bandwidth 
with reasonably steep slopes and good out of band blocking. By 
shifting the filter toward the short wavelength side of the line 
(i.e., < 8683 A) the filter bandwidth may be minimized. The 
filter will be designed for the specifics of the optical design. 
Figure 4.8 shows the wavelength shift as a function of angle and 
the approximately minimum filter bandpass attainable for the 4- 
cluster configuration. 



FIGURE 4.8 FILTER HALF POWER BANDWIDTH VERSUS 
MAXIMUM ANGLE OF INCIDENCE 
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It is evident that mean angles of incidence greater than 4 
deg will not allow the system oandpass requirements of 5 A to be 
met. The peak transmission of the 5 A filter may be as much as 
40%. The interference filter was chosen for the baseline 
configuration . 

Birefringent Filters The classic filters in this category 
are the Lyot-Ohman or Sole filters. A Sole filter was chosen for 
analysis since it is much more compact for the same FOV and 
bandwidth. 

The Sole filter can easily achieve a bandpass from 1.5 to 3 
A over incidence angles of + 10 deg. Consequently a single 
filter may be used even in a wide FOV collector optics train to 
suppress background radiation. The estimated peak transmission 
of a 5A filter for LITMAP is 25%. Figure 4.9 illustrates the 
design of a Sole filter. The use of Sole filters in a narrow FOV 
telescope is simplified since the filters are smaller and the 
angles of incidence are less. Under these conditions it would be 
relatively easy to achieve 1.5 A bandpass. 



FIGURE 4.9 SOLC FILTER 
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4.2.5 Transmission 


The optics will be antireflection coated with "V" coatings 
on each surface. These coatings will assure an inband 
transmission of 0.99 per surface. For the baseline design the 
overall transmission would be (0.99)^ - 0.3. When coupled with 
an interference filter transmission of 50%, the system 
transmission will be s 40%. In the case of the Sole filter the 
transmission will be approximately 22%. 


4.2.6 Baseline Performance Summary 


The system performance of 
Clear Aperture 
Effective Focal Length 
Pixel Size 
Resolution 
Fields-of -View (4) 

Filter Bandwdith 

Transmission 

Filter 

Configuration 


the baseline optics design is: 
9.7 cm 
9.7 cm 
320 urad 

< 140 urad 

7.4 deg x 7.4 deg 

< 5 A 

1 40 % 

Interference 
Four cluster 



4.3 DETECTOR 


Detector selection is a key element of the L1TMAP design 
process. It is here that the twin sets of constraints from the 
optical side and the electronic side converge, and it is from 
here that the logic of the design must flow, upstream towards the 
telescope and the physical appearance of the instrument and 
downstream to the signal processing and data formatting. There 
are some fairly general arguments that lead to the choice of a 
CCD as a detector, and some specific arguments that indicate the 
use of multiple CCD's to cover the full Earth disk. These are 
discussed in the following sections. 

4.3.1 Selection of Detector Type 

First, the minimum signal of interest, the worst case solar 
background, the size of a pixel on Earth and the satellite 
altitude are given. Selection of the pixel size at the detector 
fixes the focal length of the telescope. Placing a reasonable 
limit on the numerical aperture, say f / 1 , fixes the telescope 
diameter. Physical and manufacturing considerations described in 
Sections 2, 5, and Appendix A, constrain the optical filter 
bandwidth and losses. At this point in the logic the number of 
signal photons available to a detector that corresponds to the 
minimum lightning event is fixed as is the number of photons per 
second at a detector pixel due to the maximum solar background. 
For a detector pixel size of 30 um, an optical bandwidth of 5A 
and an overall optics transmission of 0.4 these numbers are 3952 
photons of signal and 1.5 1 0 8 photons/sec of background. It is 
the high magnitude of the photons per second for solar background 
that eliminates whole classes of detectors. 

For a fixed telescope f/number, changing the detector pixel 
size changes the telescope size. The number of photons received, 
both of signal and background, is proportional to the square of 
the detector pixel size. For the LITHAP application, photon 
statistics are the limiting factor In attaining a useable 
signal/noise ratio and that in particular the instrument can be 
made to work with a pixel size of 30 micrometers and not with a 
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pixel size of 15 micrometers (c.f. Appendix B). This result will 
have important consequences later. 
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Figure 4.10 shows a family tree of detector alternatives. 
"FP" on this figure stands for focal plane. The first level of 
classification of sensors divides clocked or integrating 
detectors from random event or non-integrating detectors. 
Integrating detectors, e.g. CCD's or vidicons, have a definite 
frame time: incoming photons free charge carriers which 
accumulate until the next periodic readout time. Non-integrating 
or random event detectors read out only when an event occurs. A 
one-dimensional example of this would be a photomultiplier in a 
photon counting configuration or (to sense lightning events 
instead of photon arrivals) any single detector/preamp followed 
by a high pass filter and pulse discriminator. 



• FP PROCESSING 

• FP READOUT 

OPTICAL: LOCATION 
ELECTRICAL: INTENSITY 


HIGH DATA RATE FROM FP 


LOW DATA RATE FROM FP 


FIGURE 4.10 GENERIC 0ETECT0R ALTERNATIVES 
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A true two dimensional lightning event detector with the 
background/lightning separation performed on the focal plane, is 
not now available and a significant development effort would be 
required to provide it. Much effort is currently being expended 
on the background suppression problem, particularly for infrared 
arrays, but no devices are planned with the size and frame rate 
needed for LITMAP. The following detector concepts were 
evaluated for use with off-focal-plane background suppression. 

CTD: Charge transfer devices (CTD’s on the chart) and 
vidicon-like devices are integrating devices. Silicon target 
vidicons have similar optical and detector properties to CCD’s 
but do not have the capability of multiple parallel readouts to 
keep the data rate down and are not in general suited to 
millisecond frame time applications. The only charge transfer 
device that might be considered for LITMAP besides a CCD is a 
CID, or charge injection device. A CID, with its ability to 
randomly access pixels, might be useful in a lightning mapper 
that only looked at a small, electrically steerable portion of 
the field of view, but its higher noise makes it inferior to a 
straight CCD for full time processing of the full Earth disk. 

MAMA: A two-dimensional example of a photon counter is the 
Multi-Anode MicroChannel (MAMA) system (Figure 4.11). It 
consists of a photocathode, a microchannel plate electron 
multiplier and two orthogonal sets of stripe anodes connected to 
coincidence logic. The MAMA system is an extremely effective 2D 
photon counter, but it counts photons, not lightning events and 
the maximum counting rate per pixel of less than 650 per second 
is many orders of magnitude short of the 220 million per second 
background rate. 
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FIGURE 4.11 MULTI -ANODE MICROCHANNEL ARRAY 

Anamorphic: The anamorphic concept has an appealing 

simplicity (Figure 4.12) and consists of two or possibly three 
linear arrays each with a cylindrical lens. A vertical or 

horizontal stripe from the field of view is focussed onto a 
single detector element. Each detector element has a preamp, 

high pass filter (to perform the background supression) and 
threshold circuit. The (digital) signals from the two axes are 
correlated to sense the location of a lightning strike. Because 
each detector element sees a whole stripe’s worth of solar 
background but only a pixel's worth of lightning, the 
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FIGURE 4.12 ANAMORPHIC OPTICAL SUBSYSTEM 


signal/noise ratio required for this scheme to work is quite 
high. This high signal to noise ratio could be attained in a 
near earth (e.g., Shuttle) orbit but is not possible at 
synchronous orbit. It is difficult to suppress even one pixel's 
solar background from synchronous orbit, much less 800 times 
that. 

o From the above brief discussion, the most promising 
detector for full Earth daytime lightning coverage from 
geosynchronous orbit is a CCD. These devices are now 
discussed further. 
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4 . 3 . 2 Frame Transfer Methods OF POOR QUALITY 

Fra me Transfers Figure 4.13 shows schematically how a 
typical CCD is organized and Figure 4.14 is a photograph of a 
CCD. This CCD is a frame transfer device; there is an imaging 
area which is placed at the detector focal plane and there is a 
separate image storage area. Photoelectrons accumulate in the 
imaging area during an integration time, then the entire image 
(stored electrons) is rapidly shifted intact to the image storage 



FIGURE 4,13 ORGANIZATION OF 5040 AREA ARRAY 

S> 
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FIGURE 4.14 WESTINGHOUSE 5040 

section. The image is then read out sequentially, pixel by 
pixel, a line at a time during the next integration time. Thus 
during an integration time one frame is being formed in the 
imaging section and the previous frame is being read out from the 
image storage section. A typical shift time for frame transfer 
devices is 0.3 milliseconds. This is an appreciable fraction of 
the 4 millisecond integration time needed for LITMAP and will 
cause some image smear of the stronger lightning flashes when 
they occur during the transfer interval. 

Line Transfer - There are also line transfer and interline 
transfer CCDs. Line transfer devices have no on-chip memory and 
must use an external shutter to block the image while they are 
being read out. These are useful only in astronomical 
applications and the like where the significant readout time can 
be tolerated. Interline transfer devices have columns of opaque 
memory cells interleaved with columns of photosensitive cells. 
On command the accumulated electrons are shifted sideways one 
line into the interline memory columns where they are read out 
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during the next integration time. This type of CCD has much to 
recommend it for high frame rate applications: for frame 
transfer, the image moves over only one pixel instead of down a 
whole column and so transfers faster in the ratio of the number 
of pixels per column. The disadvantage of interline transfer 
devices is the lower effective quantum efficiency (0.12) 
resulting from having half the sensing area opaque for the frame 
storage function. Recently, rows of tiny cylindrical lenses 
have been manufactured with the sensor chip to redirect most of 
the focal plane light onto the optically sensitive portions 
resulting in increasing the effective detection quantum 
efficiency. 

4.3*3 Focal Plane Configuration 

Four CCD focal plane configurations were addressed to 
provide an assessment of alternate approaches for satisfying the 
LITMAP objective of obtaining lightning data over a full disc 
from geosynchronous orbit. These configurations ranged from 
severely pushing the existing technology to an approach which 
utilizes existing, unmodified hardware. Approach B is the 
baseline used for the concept definition and costing tasks of the 
study. Briefly these configurations are: 

A: Two Chip - 800 X 800 Array - Modified 

This is the most discussed approach using two 
sensor arrays such as the TI 800 X 800 array that is 
used by astronomers and by NASA for the Galileo 
project. A single simple prism split could be used to 
avoid any loss of information along the mating edge. 
No array exists which can handle the data rate 
required. Primary modification of existing devices 
would be array partitioning and the adding of output 
buffers and amplifiers to provide between 8 to 16 
outputs per array, each operating in the 10 MHz regime 
to provide a two to four millisecond frame readout 
time. System analysis indicates that a pixel size 
considerably greater than the 15 microns of the 
existing 800 X 800 TI device is needed to handle the 
daytime background and to provide the dynamic range. 
TI has indicated that the yield for such a large device 
would be so low that its production could not be 
justified. 

Mosaic: An alternative for the single focal 
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plane is to further sub-divide the focal plane and butt 
the arrays together into a single mosaic. Loss of 
significant information and the problem of getting a 
narrow band filter for the single wide FOV are 
significant added complications. 

Itek/Bell Northern has been able to precisly cut 
and butt mount (along two opposite edges) sensor arrays 
to form a long continuous array. This technique has 
been extended to three sided buttability by several 
sensor manufacturers with RCA claiming that they should 
be able to do four sided butting with their PtSi sensor 
technology. In all cases to date, the dead space is at 
least 50 microns. For LITMAP this results in loosing 
one to three pixels. The affect on a mosaic sensed 
image is the loss of information along strips which are 
between 10 to 30km wide at the sub-satillite point. 
Depending on the mosaic arrangement and satellite 
location, this loss of information could be 
intolerable. An operational approach could be used 
which slightly repositions the sensor FOV whenever 
coverage under the normally missed areas is desired. 

Although the mosaic approach could be used to 
cover the full FOV with a single optical system with no 
auxiliary image splitting, there remains the problem of 
providing an optimum narrow band filter to cover the 
wide FOV as discussed eleswhere. The mosaic is a 
viable alternative and must be considered during 
detailed instrument design. 

B: Four Chips - TI, RCA, Fairchild - Modified 

Optically subdividing the total FOV into four 
sections relaxes the FOV and pixel cell size 
constraints to the point of being cost effective and 
compatible with the present optical and sensor 
technologies. Additional array partitioning is still 
required to accomodate the necessary data rate. Such 
partitioning, although fairly straight forward, will 
still require on the order of a million dollar 
development activity with about a year lead time. This 
is the baseline approach that is primarily focused upon 
in this final report. 

The FOV splitting can be done with a single 
telescope coupled with four way prism image splitting 
or by using four independent telescopes for the four 
sub arrays. The required low F/No makes the prism 
splitting very difficult hence the four telescope 
approach is baselined. (A mosaic and single focal 
plane was discussed in the first configuration above.) 
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C: . Four Chips - RCA - Unmodified 

A RCA chip exists which can provide the required 
exposure time via the following procedure. The basic 
frame time through the existing output port is 17 
milliseconds. The image section can be exposed for 13 
milliseconds and then dumped to a drain at the opposite 
edge away from the output mux. The next 4 millisecond 
exposure is transferred to the on chip frame storage 
for output during the next 17 milliseconds. The 
dumping/output cycle can be repeated giving a 3 
millisecond (or any other desired integration time) 
sample every 17 milliseconds. 

This time sampling approach can provide statistics 
on total energy/stroke distribution but there will be 
no means of separating first return stroke data from 
other data. For an operational system such as a severe 
storm indicator, this mode of operation may be 
satisfactory. For scientific objectives, the relative 
roles of the strokes (first versus subsequent) will not 
be retrievable. 

A potential variation of this time sampling 
approach would be to, during the readout cycle, rapidly 
transfer m rows (One or two MHz rate) to the output 
mux which will result in a loss of this data. The next 
n rows are transfered and processed at the normal rate 
(sub MHz row transfer). The remaining rows would be 
rapidly dumped, again with the total cycle time being 
the allowable 4 milliseconds (or whatever the selected 
integration time). The same n rows could be processed 
for a period of time (a few seconds would give hundreds 
of sampling cycles) then shifted to the next n rows. 
Repeat measurements over the same area would occur 
every few minutes providing a time sequence for each 
active area hence providing full coverage on a time- 
share basis. This approach would preserve the 
distinction between first and subsequent strokes which 
has scientific value. 

D: Single (or two) Small Arrays With Mechanical Scan - 
Unmodified 

A Westinghouse 2 X 100 x 100 chip exists which can 
provide the frame time required. Since only a few 
percent of the total field of view will contain active 
storms at any time, a smaller array combined with a 
mechanical scan mechanism can be used to cover the 
desired area. A given sub area can be observed for 
several seconds giving a few hundred frames and then 
another area is measured. On the order of 10C 
positions would cover the full field of view, hence a 
recycle over a given area can be done every few 
minutes. Alternately, a given storm can be fixed upon 
and followed over a period of time with only an 
occasional scan over the full area to determine if 
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other severe storm areas are developing. 

The above summarized four focal plane concepts which can provide 
lightning data. Concept B will provide the best full time 
coverage at the expense of special sensor development. Approach 
D provides an off-the-shelf concept which can provide good storm 
coverage at the expense of an added mechanical tracking 
mechanism. Approach B was chosen for the baseline technology 
assessment . 

o Approach B was selected for the concept definition and 
costing tasks of this study. 

4.3.4 Required Modifications 

Table 4.1 illustrates a signal/noise calculation for a RCA 
CCD type SID 501D. The three left-most columns of the day and 
night tabulations on the bottom half of the page provide the 
minimum detectable and maximum non-saturating lightning intensity 

tabulated as a function of frame time. The units of the 
lightning intensity are uJ/m2_sr. These numbers can be compared 

with the MSFC NASA/U2 data taken in the summer of 1982: 

Lightning Intensity Percent of Events 

(Max integrated Below this 

radian?e/event) Intensity 

uJ/m^ sr 


1.8 

1.6 

4.2 

7.2 

5.6 

12.0 

7.5 

22.3 

13.3 

50.9 

42.2 

80.2 

56.2 

89.7 

100.0 

98.4 


A look at Table 4.1 shows that a frame time of two 
milliseconds gives excellent performance at night, losing only 
about 2 percent of lightning events at the low end and none at 
the high end. The daytime performance for the same integration 
time is adequate: a loss of 12 percent of lightning events at 
the low end and again, none at the high end. 

The picture with a four millisecond frame time is a little 
less optimistic: the system performance is essentially unchanged 
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TABLE 4,1 RCA SIO 501D Performance 


* * * lxtmap s/n * * ♦ 


This is a calculation of tha LITMAP minimum detectable signal 
and aaNiaua non-saturating signal expressed in terms of lightning 
intensity at tha cloud tops in (uJ/m : ir), 


Required S/N ratio « 3.9 

Pixel size on Earth (edge of square) (km) * 10 

Effective source area (km'2> *• 100 

Average stroke duration (ms FWHH) * .34 

Satellite altitude (km) * 33600 

CCD type - RCA SZD 30 1 D 

Pixel size of CCD (urn) * 30 

Quantum efficiency * .28 

Readout noise (rms els) * 150 

Full Mali (electrons) * 230000 

Telescope f /number * 1 

Telescope focal length (cm) * 10.66 

Telescope diameter (cm) » 10.68 

Solar irradiance <W/m*2 sr) * 301 

Cloud albedo • .8 

Optics transmission, including filter - .5 

Wavelength (nm) » 868.3 

Filter bandwidth (nm) ■ .6 


Day (worst case sun ) : 


Frame 

Lightning 

Pulse 

Solar 

Hin 

Bkgrnd 

Readout 

Total 

Time 

Hin 

flax 

Splitting 

8k grnd 

Signal 

Noise 

Noise 

Noi se 

(ms) 

(uJ/m~ 

2 sr) 

Loss 

(els) 

(els) 

(res 

. electrons) 

0.23 

6.9 

1385 

0.39 

13614 

1169 

123 

130 

198 

0.30 

3.0 

791 

0.64 

31227 

1383 

177 

150 

233 

1.00 

4.9 

329 

0.82 

62435 

1737 

230 

130 

294 

2.00 

3.8 

318 

0.91 

124909 

2283 

333 

130 

387 

4.00 

7.3 

0 

0.96 

249818 

3096 

300 

130 

323 

8.00 

10. 1 

« 

0.98 

499636 

4281 

707 

130 

726 

16.00 

14.0 

• 

0.99 

999273 

3981 

1000 

130 

1014 

* Solar background 

1 exceeds 

ful 1 wel 1 

capacity 

* a 




Night: 

Frame 

Lightning 

Pul se 

Solar 

Hin 

8k grnd 

Readout 

Total 

Timm 

Hin 

Hax 

Splitting 

Bkgrnd 

Si gnal 

Noi se 

Noise 

Noise 

(ms) 

(uJ/m‘ 

*2 sr) 

Loss 

(els) 

(els) 

(rms electrons) 

0.23 

3.3 

1477 

0.39 

0 

903 

0 

130 

153 

0.30 

3.3 

903 

0.64 

0 

903 

0 

130 

153 

1.00 

2.3 

703 

0.82 

0 

903 

0 

130 

133 

2.00 

2.3 

636 

0.91 

0 

903 

0 

130 

133 

4.00 

2.2 

606 

0.96 

0 

903 

0 

130 

133 

8.00 

2.1 

392 

0.98 

0 

903 

0 

130 

133 

16.00 

2.1 

383 

0.99 

0 

903 

0 

130 

133 
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at night, but in the daytime the CCD is saturated by the solar 
background. Even if the full well capacity were higher, the low 
end threshold has risen to the point of losing 20 percent of the 
lightning events. Hence the framing time needs to be kept as low 
as possible and the full well capacity as high as possible (see 
Appendix A for further analysis). 

Commercial CCD's have frame rates appropriate to TV 
applications or slower, i.e. 50 or 60 per second or less, and can 
be driven less than a factor of two faster than that before 
significant performance degradation occurs. A frame rate of 250 
frames per second (4 millisecond integration time) is required. 
The restriction is the speed at which a CCD can be read out which 
is limited primarily by the output amplifier. Image partitioning 
with the data sent to multiple output ports resolves this 
bottleneck. The manufacturers of the sensors have proposed 
various schemes for doing the partitioning (c.f. Technology 
Assessment, Section 5). 

The electronics to deal with the multiple data streams 
emanating from the detector must be made in multiple parallel 
versions also. This electronics and how the background 
subtraction is to be implemented are discussed in the next 
section. 

Table 4.2 is a signal/noise calculation for the device 
originally proposed by the Optical Science Center, University of 
Arizona, for a lightning mapper, the Texas Instruments 800 x 800 
array. The threshold levels are above 25 uJ/m 2 - sr due to the 
5.4 cm aperture which is determined by the small 15 micron pixel 
size. Comparison of this 15 micron device results with the 30 
micron device results above (c.f. Table 4.1) indicates that a 
large pixel size is required to meet the LITHAP requirements 
(c.f. Appendices A, B, and Sections 2 and 3 for further details). 

o Integration time must be kept low (4 milllsec nominal). 

o Array partitioning is required to keep the port 
operating frequency within technology limits. 

o Large pixel sizes are needed to attain low threshold 
capability. 
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TABLE 4.2 TI 800 x 800 Performance 


* * * LITHAP S/N * * * 


This is « cftl cul fttlon of the LITHAP Minimum detectable signal 
ftnd maximum non-saturating ftign«l exprassed in terms of lightning 
i ntftifti ty ftt thft cloud top« in (uJ/m~2 sr). 

Required B/N rfttio - 3.9 

Pixel ftixft on Earth (edge of square) (km) ■ iO 

Effective ftourcft ftrftft <km A 2> » 100 

Average stroke duretion <m« FWHH) • .34 

Satellite altitude (km) ■ 39400 

CCD type • Tex aft Instruments SCO x SOO 

Pixel size of CCD (urn) - 13.2 

Quantum efficiency ■ .13 

Readout noise <ro« els) * 130 

Full well (electrons) ■ 109000 

Telescope f /number » i 

Telescope focal length (cm) « 5.4112 

Telescope diameter (cm) • 3.4112 

Solar irradiance (H/m*2 s r) * 301 

Cloud albedo ■ .8 

Optics transmission, including filter • .3 

Wavelength Cnm) * 866.3 

Filter bandwidth <nm) » .6 


Day (worst case sun): 


Frame 

Li ghtmng 

Pul se 

Solar 

Hin 

Bkgrnd 

Readout 

Total 

Time 

Hi n 

Hax 

Splitting 

Skgrnd 

Signal 

Noise 

Noise 

Not se 

(ms) 

(uJ/m' 

‘2 Sr) 

Loss 

(els) 

(els) 

(rms electrons) 

0. *?3 

46.3 

3311 

0.39 

1S41 

939 

43 

130 

139 

0.50 

29.3 

3192 

0.64 

3722 

973 

41 

130 

163 

1.00 

24.4 

2404 

0. 82 

7444 

1039 

94 

130 

176 

2.00 

24.7 

2007 

0.91 

14989 

1199 

122 

150 

196 

4.00 

27. • 

1410 

0.94 

29773 

1344 

173 

130 

232 

8.00 

33.9 

9S2 

0.99 

99930 

1709 

244 

130 

2S9 

16.00 

43.9 

a 

0.99 

119101 

2239 

343 

130 

379 

• Solar background exceeds full well 

capacity. 



Night: 









Frame 

Liohtning 

Pulse 

Solar 

Hin 

•kgrnd 

Readout 

Total 

Time 

Hin 

Hex 

Splitting 

Skgrnd 

Si gnal 

Nol«* 

Noise 

Noise 

(ms) 

(uJ/m~2 er> 

Loss 

(els) 

(els) 

(rms electrons) 

0.29 

44.7 

9404 

0.39 

0 

903 

0 

190 

193 

0.90 

27.4 

3309 

0.44 

0 

903 

0 

190 

193 

1.00 

21.4 

2990 

0.92 

0 

903 

0 

190 

133 

2*00 

19.3 

2329 

0.91 

0 

903 

0 

190 

13 T 

4.00 

IS. 3 

2219 

0.94 

0 

903 

0 

190 

13L 

9.00 

17.9 

2144 

0.99 

0 

903 

0 

190 

133 

14.00 

17.7 

2140 

0.99 

0 

903 

0 

190 

193 
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4.4 SIGNAL PROCESSING 


The signal content of each pixel consists of various 
components of which .<• desirable component is the lightning 
stroke intensity integrated over the sample time. The unwanted 
components consist of integrated background irradiance and the 
detector noise. Etch pixel signal arriving from the detector 
must be tested for the presence of a lightning stroke. The 
detection, extraction and measurement ~>f the lightning signal 
from the detector output data stream is the primary function of 
the data processor. The data processor will perform the 
following functions: 

a. Detect lightning stroke with time and location tag 

b. Measure lightning energy 

c. Format and return the data 

The baselined 800 X 800 CCD imaging portion of t'le array 
consists of 6.4 X 1 0 5 pixels which must be individually 
integrated for the presence of lightning at the sampling rate of 
4 milliseconds or less per frame. This transforms to a total 
data rate of at least 160 mega samples per second for the total 
detector array. This data rate through a simple port exceeds the 
present electronic processing as well as the CCD technologies. 
The data rate can be reduced to an acceptable level for 
processing by using multiple output ports from the detector array 
(see previous section). Table 4.3 gives the various combinations 
of array size, output port number, sampling time and the 
resulting sample rate. Each output port will then be serviced by 
a signal processor circuit. The extracted event data from each 
data processor, however, can be streamlined into a single data 
formatting circuit with reduced data handling rate capability 
because the lightning occurance rate is sufficiently low. This 
parallel processing approach reduces the complexity of the system 
without sacrificing the performance. 

4.4.1 Design Requirements 

The data processor must be designed to perform the three 
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FIGURE 4.15 LITMAP SYSTEM CONCEPTUAL BLOCK DIAGRAM 


primary functions listed in the previous section. In order to 
facilitate the evaluation of various design approaches and 
concepts, baseline design requirements were established. The 
baseline system block diagram in Figure 4.15 shows the primary 
functional blocks for a data processor servicing one of the 
output ports of the CCD. The total size and power requirement of 
the system is very much dependent on the number of output ports 
which results in the duplication of this data processer section. 
Based on past experience each data processor is estimated to 
occupy approximate half of a 7" X 7” circuit board with about 
eighty 16-pin dual-in-line packages. Therefore a system with 16 
output ports can be fitted into less than one cubic foot volume 
which is comparable to other flight systems of similar 
complexity. From Table 4.3, a 16 channel system with an 800 X 
800 CCD device at 4 ms sampling interval is equivalent to an 
effective data rate of 10 mega samples per second. A reasonable 
power budget based in similar flight hardware, is in the range of 


47 




ORIGINAL page is 
OF POOR QUALITY 


♦) 


TABLE 4.3 Array Size and Data Rate Trade. 


MMT SIZE 

vmi interval 
(NSEC* 

NO. OUTPUT 
PORTS 

SAMPLE DATA RATE 
PER PORT (SPS) 

BACKGROUND PEMORT 
SIZE PER PORT 
(SAMPLES) 

PROCESSING TIME 
PER SAMPLE 

(NSEC) 

40 x ioo 

3 

16 

13.3 I 10 6 

40000 

75 

m x «oo 

3 

32 

6.7 X 10* 

20000 

ISO 

m x no 

4 

14 

10 X 10 6 

40000 

10 

40 1 100 

4 

32 

5 X 10 6 

20000 

200 

40 IM 

5 

16 

8 X 10 6 

40000 

121 

00 K 80 
• 

5 

32 

4 X 10* 

20000 

250 

00 X 00 

3 

16 

10 X 10* 

30000 

10 

00 x on 

3 

32 

5 X 10 6 

15000 

20 



00 X 80 5 


» ( X 10* 

3000 

167 


80 X 80 5 


r j x io* 

1500 

333 


100 X 10 3 


1 3.3 X to* 

1000 

30 

j 

r 

10 X 10 3 


« 0.8 X 10* 

250 

120 


10 X 10 5 


1 2 X 10* 

1000 

50 

* 

10 X 10 5 


« 0.5 X 10* 

2500 

200 

4 ' 

0 . 

T 

30 to 50 watts. 





* 

In summary, 

the 

baseline electronic 

i subsystem 

requirements 

- 

are; 





a 

1 . Data rate 

- 

10 mega samples 

per second 

per output port 


2 . Power 

- 

30 to 50 watts 




3 . Size 

- 

approximately 1 

cubic feet 


t 

4 . Weight 

- 

approximately 60 

pounds 



4 . 4.2 Background Subtraction 

The lightning mapper system is required to detect lightning 
flashes during day and night. Since the background conditions 
vary drastically between day and night and over cloud top, land 
and sea, it is necessary to design the system to handle the worst 
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case background illumination conditions. 

There are essentially three methods for extracting weak 
signals from a strong background. These are spectral, temporal, 
and spatial filtering. The spectral and temporal filtering may 
be implemented at the optical and analog signal stages of the 
process. They may be considered predetection processing 
techniques. Spectral filtering is essential to the performance 
of the system when the clouds are sunlit. The use of an optical 
filter has been discussed in detail in an earlier section and 
will not be presented here. 

Spatial filtering uses an averaged background level from 
some adjacent pixels which is subtracted from the incoming signal 
in order to isolate possible lightning signals. This technique 
requires that the background level be relatively uniform over 
adjacent areas which often is not a valid assumption. Therefore 
the use of spatial filtering for background subtraction will not 
be considered for the LITMAP application. 

Temporal filtering depends on the background remaining 
static for two or more successive frames. Successive frames may 
be subtracted to elliminate the background signal from each pixel 
with the drawback that the overall noise will increase by a 
factor of s/2. A running average technique for the background 
statistically can reduce the s/2 factor significantly. Further 
discussion of this technique will be presented in a later 
section. 

4.4.2. 1 Analog vs Digital 

T he baseline conceptual block diagram of the lightning 
mapi. > system is shown in Figure 4.15. Temporal background 
subtraction is employed. The present frame pixel signal S(i) is 
being fed into the background subtraction amplifier. The 
subtrahend B(i) is equivalent to SC i— 1 > (S(i) from the previous 
frame) plus some contribution from the past due to averaging. 
This scheme can be realized by using analog or digital hardware, 
or a hybrid. These approaches will be analyzed against the 
technology and design requirements and then, a data processor 
system design concept will be discussed. 
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4. 4. 2. 2 Digital Approach 

The digital background subtraction system can be visualized 
from Figure 4.15 by installing the A/D converter for the digital 
system as shown. The input signals from the focal plane CCD 
detector are preprocessed by the A/D to provide digital data. 
Subsequent data processing such as multiplication and background 
subtraction is carried out in digital circuitry. The three 
important circuit components to be considered in this digital 
approach are: the A/D converter, the background memory and the 
associated logic circuitry that controls the processing. In 
order to develop a feasible flight system the critical factors 
that need to be considered are size, weight, power and technology 
availability. 

The A/D converter is a very important component in this 
approach. Based on previous a.ta and calculations, the 
background irradiance can occupy as much as 20% of the i D full 
well capability. IF the minimum lightning signal is taken into 
account, the A/D converter will require better than nine bits of 
resolution and the background memory will require eight bits for 
background information storage. The A/D resolution and range 
requirement can only be satified simultaneously by flash type A/D 
converters which have rather high power requirements (on the 
order of one watt). Lower power devices are available with less 
bit resolution but would sacrifice resolution and therefore can 
only be considered as an alternative. Another way to reduce the 
A/D resolution is by front-end signal compression using a 
logarithmic amplifier. Application of the log amplifier at this 
point, however, will induce large inaccuracy in the system where 
background signal is a significant portion of the entire useful 
signal range. This inaccuracy is worse at low level lightning 
signals and consequently will tend to reduce the quantity and 
quality of useful data. The signal compression technique, 
however, will be useful after the background subtraction to 
reduce the required resolution of the A/D converter for the 
analog approach that is discussed in a following section. 

The data processing technique that is proposed in this 
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system consists of background averaging, limiting, subtraction 
and adaptive threshold adjustment. In order to realize these 
functions digitally, high speed digital processing circuitry and 
components have to be utilized which tends to significantly 
increase the power consumption of c,r = system. 

The background storage memory is a digital storage device 
and can be a module made up of an array of random access memory 
devices (RAM) with the associated logic control circuitry. The 
proposed system requires that the signals derived from the pixels 
within one frame be processed within the next frame of time. 
This results in less than 100 nsec data time per pixel based on 
10 MSPS data rate. In order to store the current digitized frame 
information, the background memory in a digital system will 
require a storage capacity of 40000 X 8 bits (for 8 bit A/D 
conversion) at an access time of much better than 100 nsec. 
Commercially available bipolar or highspeed M0S devices generally 
can operate with better than nsec access time for read or write 
cycles. A module of 40000 X 8 bit memory will typically consists 
of 20 4k X 4 bit RAM chips in addition to the associated logic 
circuit for addressing control, and signal processing. When 
translated into 16 parallel output channels the total electronics 
size, weight and power will exceed the proposed design 
requirement by a wide margin. Therefore, alternate methods or 
devices have to be sought. 

One method of reducing the per channel power requirement is 
to use low power high density M0S type RAM memories. Static RAMS 
with 8k X 8 memory size and 1/4 watt per chip power dissipation 
are available but are generally slower. The most straight 
forward approach to oversome the slow speed is to increase the 
number of output ports along with the number of data processors. 
If the 10 MSPS data rate is reduced to 5 MSPS by increasing the 
number of output ports to 32 for a 800 X 800 CCD array, the 
processing time is increased to 200NS and the background memory 
size is reduced to 20000 elements per data processor. This 
translates into a total of 3 memory chips and 3/4 watts power 
requirement per data processor. The power reduction however, is 
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offset by an additional 6 data processor modules. An alternate 
wayof utilizing these lower power, higher density devices is to 
multiplex the data to two or more memory modules per channel to 
increase the effective storage area at a slower speed. However, 
all of this will require additional logic control circuitry and 
will increase overall size and power. 

To summarize, the digital approach is best if size and power 
are not constraints. Since the objective is to design a system 
possessing reasonable size and weight, digital method appears to 
be a questionable approach. The future trend is toward high 
speed, lower power and high density. It is possible that in the 
future a system based on a digital approach would be more 
compatible with the power constraints wnile providing the 
necessary processing speed. 

4. 4. 2. 3 Analog Approach 

An analog approach to background subtraction can be realized 
by storing the previous frame pixel signals in a CCD memory. A 
CCD device developed for the focal plane can be used for this 
purpose by masking out the photo sensitive area. On chip input 
and output ports equal in number to the focal plane output ports 
can be implemented. In addition the approach is a cost effective 
way of utilizing existing technology. Analog approach involves 
preprocessing the time multiplexed analog signal with fast analog 
circuitry. This type of circuitry currently exists and 
therefore the approach represents a lower risk than the digital 
approach. The total size, weight and power requirements are also 
substantially lower than for the digital approach. Further 
discussion of this technique will be presented below as part of 
the baseline approach. 

4.4.3 Baseline System 

Figure 4.16 shows the proposed data processor system block 
diagram. The focal plane sensor array consists of one or several 
large array CCD detectors, each with built in multiple-output 
ports. The system will be operated in the staring mode with 
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FIGURE 4.16 PROPOSED LITMAP ELECTRONICS SYSTEM BLOCK DIAGRAM 


frame rate in the two to five milliseconds duration to optimize 
SNR. Detail discussion of the device and the SNR were presented 
in earlier sections and will not be repeated here. Each output 
port is serviced by its own analog data processer (DP). 
Background subtraction and threshold detection are carried out 
within the DP. The subsequent lightning intensity data, along 
with the time and location tags, will be fed to a single data 
formatter (DAF). The extremely low lightning occurance rate 
allows the use of a single DAF approach without compromising 
system throughput while minimizing system complexity and 
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maximizes system efficiency. This multiple output-port approach 
also permits scaling the system (number of output ports) to match 
the array size with minimal modification of the baseline system. 
Example of this would be a small array focal plane lightning 
mapper system with a single DP that could be flown on Shuttle 
mission in lower earth orbit. A detailed description of the 
system is presented in the following sections. The resulting key 
features for the signal processing are: 

1. Multi-port detector array 

2. Multi-channel analog data processor 

3. Background subtraction by temporal technique 

4. Background smoothing by running average technique. 

5. Post detection digitizing using A/D converter 

6. Single data formatting circuit. 

4.4.3. 1 Data Processor 

Outpu t Buff er. The proposed data processor (DP) as 
described earlier is based on an analog approach. The CCD output 
signal is buffered and the data sampled and held. This S/H 
function is utilized to sample the data at the appropriate time 
to avoid transients and switching glitches. The buffer also 
serves to extract the signal from the DC bias voltage of the CCD 
output amplifier. Commercially available instrumentation and 
operational amplifiers and sample and hold amplifiers can be used 
to implement this curcuit. However, at 100 ns processing time 
per pixel, the bandwidth requirement for these devices are well 
above 10 MHz. Commercially available devices which can operate 
to this specification generally have high power requirements 
relative to the total system power budget established in Section 
4.4.1. Consequently, it is necessary to design circuitry based 
on discrete components or hybrid elements which also provides 
additional size savings. 

The buffered signal is fed into the background averager and 
the background subtraction amplifier for further signal 
processing. The description and design approach of these 
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functional blocks will be discussed next. 

Background Subtraction . The function of the background 
subtraction circuitry is to isolate lightning signals from the 
background and thus reduce the dynamic range required for the A/D 
converter. Averaging techniques are used to increase system 
responsitivity and performance. A simplified functional block 
diagram is shown in Figure 4 . 17 . Subtraction is carried out in 
the subtraction amplifier with inputs S(i) and B(i). S(i) is the 
buffered output from the detector and B(i) is a conditioned 
i derivative of S(i) from previous frames. 

i The input signal S(i) is fed to the amplitude limiter which 

i prevents excessive lightning signal amplitude passing to and 

contaminating the background data. The limiter provides a coarse 
limiting of the signal pulse with further smoothing provided by 
time averaging the background samples. An adaptive feedback 
based on the average background is used to optimize for both day 
and night operations. 

The background averager consists of the multipliers k and 
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(1-k), the adder and the background memory. The background 
memory is essentially the same type of CCD device used in the 
focal plane detector. The function of this device is to store 
the analog data arriving from the adder, and to output on a 
first-in first-out (FIFO) basis the stored data for the 
background subtraction and for the averaging routine. The 
coefficient k is confined between zero and one. In general, the 
sum of the two multipliers k and (1-k) should be less than unity 
to prevent the system from saturating. The input to the 
background memory A (i) is 

A (i) = K S (i) = (1-K) B (i) 

with B (i) = A ( i — 1 ) 

where i and i-1 refer to present and previous frames 
respecitvely . 

If K=0, there is a straight throughput with no subtraction. 
On the other hand, k=1 corresponds to a simple f rame-to-f rame 
differencing. For 0 < k < 1, the subtrahend S ( i — 1 ) is a "fading 
moving average". Frames are filtered with a time constant 
inversely proportional to k, i.e., the relative contribution of 
each sample decreases with age. Let S(i) contain some noise Np* 
Then R(i) will contain some noise Nr also given by: 

Nr/Nf s (2/(2-k)) 1/2 

Some representative values of N R /Np are: 


K 

Nr/NF 

N 


0 

1.000 

•0 

No subtraction 

0.001 

1.000 

1000 


0. 125 

1.033 

7.49 

Selected N 

0.25 

1.069 

3.46 

Another ADD for 
(1-k) = 3/4 

0.5 

1.155 

1.44 

Additional ADD 

0.75 

1.265 

< 1 


1 

1.1114 

< 1 

FIFO simple two 
frame subtraction 
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where N is the effective number of frames being averaged as 
defined by 

(1-k)N = e -1 or N f 1/ In ( 1 / ( 1 -k) ) . 

Conversely 

k = 1 - e-1/N. 

As the number of background frames being averaged increases the 
noise approaches that of the original frame (assuming incoherent 
noise). Simple frame to frame subtraction (k=1) will increase 
noise by J~2, hence increasing the threshold setting. Background 
averaging (five to eight frames) permits the background to be 
removed while minimizing signal degradation. 

The overall background substraction circuitry can be 
realized in hardware by an integrated rather than a modular 
design. The multiplexers k, k-1 and adder can be implemented by 
a single summing OP AMP thus reducing the component count and 
therefore power requirements. Further size reduction can be 
achieved by hybridization of the entire background subtraction 
circuitry (excluding the background memory). It was mentioned 
earlier that the wide bandwidth requirement tends to drive up the 
power consumption of an OP AMP. It is therefore necessary to 
select the device with tl}e best power to bandwidth trade but 
still be able to meet the performance specification. Commercial 
OP AMP components will be selected where available while 
transistor circuits using single supply source will be designed 
when necessary to minimize power consumption. 

Another important design consideration is the circuit 
timing. The high data, rate requirement provides very little 
excess processing time to allow for delay, transients and 
settling. Therefore, it is necessary to minimize these timing 
errors by optimally applying high speed circuit design techniques 
such as minimizing load capacitance, path and lead length, and 
transistor saturation. 
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Background Memory The background analog memory consists of 
several CCD devices similar to those used in the focal plane. 
The devices are opaque and have electrical input ports to accept 
the analog data. The background memory and supportive circuitry 
can be built on a single circuit card. 

Input and output buffers are required for signal 
conditioning and signal level matching. On the input side of the 
memory a sample and hold (S/H) amplifier samples the output from 
the adder during the allowable time window to avoid transients. 
Amplification raises the background signal to a level which 
efficiently fills the full well capacity of the CCD memory to 
minimize the noise contribution from the analog memory. The 
output ports of the background memory have similar S/H 
amplifiers. 

Post-Subtraction Processing. Post-subtraction processing 
consists of digitizing the lightning signal and buffering it with 
the appropriate location tag for distribution to the data 
formatter. Data compression using logarithmic amplifier may be 
implemented here to reduce the resolution requirement of the A/D 
converter. A threshold detector is utilized to establish the 
false alarm rate prior to digitization thus minimizing the load 
on the A/D. The functional block diagram of the post subtraction 
processing is shown in Fig. 4.16. 

The primary function for an automatic background adaptive 
threshold control is to vary the threshold for day and night 
conditions, which occurs fairly slowly. The bandwidth 
requirement of this circuit is estimated to be much less than the 
focal plane data rate. Therefore, the design of this circuitry 
does not present any critical problems. Manual ground control 
circuitry can be implemented by feeding a digitized telemetry 
signal through a D/A converter and applying the analog voltage to 
the threshold detector circuitry. 

The output of the background subtraction amplifier consists 
of pulses coinciding with the occurrence of lightning strokes. 
The average occurrence rate is extremely low. The maximum rate 
anticipated is in the order of several strokes per frame (entire 
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focal plane). The temporal spreading of a single stroke across 
several integration intervals is treated as two lightning events 
occuring in different sample frames and, therefore, poses no 
difficulty to the post-subtraction processing circuitry. The 
spatial straddling, however, will require the post-subtraction 
processing circuitry to be able to handle bursts of data packe t 
associated with consecutive pixels. At 10 MSPS data rate the 
post-subtraction processing circuitry is required to complete 
each single event processing within 100 nsec. Hence, the 
bandwidth requirement on this circuitry is similar to the 
background subtraction circuitry except that it will be in a 
quiescent state most of the time due to the low average event 
rate. The low data rate also results in a lower powe r level than 
the similar circuitry used for the background subtraction. 

Data Conversion . The output of the threshold detector is 
a one/zero signal to indicate the presence or absence of a 
lightning signal. This signal will be used to command the A/D 
converter to initiate conversion upon positive detection of ? 
lightning signal, and to inform the control logic for timing and 
location tagging functions. 

As mentioned earlier, the detected lightning stroke 
intensity is compressed by a logarithmic amplifier and a sample 
and hold amplifier freezes the output signal from the background 
subtraction amplifier. The bandwidth requirement of this circuit 
will still be dictated by the high event burst data rate as 
discussed in the last section, hence it will be Implemented with 
operational amplifier technique using discrete components. 

The A/D converter performs the final digitizing function on 
the log-compressed signal. The dynamic range requirement of the 
signal is 100:1. It is necessary to have adequate resolution in 
order to extract useful statistical data for the lower lightning 
intensities. The log amplifier compresses the Intensity data 
without sacrificing the low level signal resolution. The lowest 
bit resolution of the A/D converter will depend on the log 
amplifier output level. In order to utilize the bit counts of 
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the A/D converter over the dynamic range efficiently, it is 
necessary to set a minimum bit resolution level. The minimum 
signal level based on a 4 msec frame time and 5A filter bandwidth 
is about 1800 electrons. Based on minimum bit resolution equal 
to 600 electrons and full well capacity of 250k electrons, a 6- 
bit A/D converter is capable of delivering an adequate 
performance. The table below lists the bit resolution for 
various compressed signal levels from 6, 7, and 8-bit A/D 

converters : 


SIGNAL LEVEL RESOLUTION AT STATED SIGNAL LEVEL 








7~-BIT 

" 8-BIT 

250k 

el 

(Full 

well) 

22828 el 

11597 el 

5845 el 

110k 

el 



10044 el 

5103 el 

2572 el 

10k 

el 



913 el 

464 el 

234 el 

1830 

el 

(Min. 

signal) 

167 el 

85 

43 el 


Each additional bit doubles the resolution capability of the 
output . 

In order to have adequate slew rate capability, the A/D 
converter must have a sampling frequency double that of the data 
rate. Presently, off-the-shelf 6-BIT A/D converters capable of 
delivering this performance are available. RCA and Siemens ere 
delivering devices requiring less than 1/2 watt. These devices 
(and their next generation spin-offs) will be the primary 
candidates selected to meet the design requirement. 

Data Buffer . The function of the data buffer (Fig. 4.18) 
is to temporarily store the A/D output data and tne location 
information from the control logic. This information will be 
extracted by the data formatter for downlink to the ground 
station. 

As discussed earlier, each data processor (channel) will 
Interface with some 40,000 pixels during each frame of sampling 
time. The location of these 40,000 pixels can be labelled by a 
16-bit word. The proposed system will cycle throug.: the 40,000 
location labels in synchronization with the sequence of signal 
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FIGURE 4.18 DATA BUFFER FUNCTIONAL BLOCK DIAGRAM 


processing. When a lightning signal is detected by the threshold 
detector, its location label will be latched into the data buffer 
memory in 8-bit words, followed by the A/D output word. Thus, 
there will be three 8-bit data words to be stored in the data 
buffer for each pixel containing a lightning signal (Fig. 4.19). 
Based on the occurance rate of lightning it is estimated that no 
more than four events will occur within the same 40,000 pixels in 
the same frame. Therefore, the size of the memory required is 
about 12x8 bit. The speed of the memory depends on the way the 
data words are latched into the memory. An 8-bit parallel data 
configuration will require an access time better than one third 
of the processing time, or better than 33 nsec. One method to 
implement this is to use two low power, high speed CMOS (such as 
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FIGURE 4.19 DATA BUFFER MEMORY STAGE FORMAT 


QMOS) FIFO memories. Due to the low average occurrence rate, the 
FIFO will essentially be at low standby power. The data 
formatter will ha die the channel tag and the time tag. Further 
discussion of the data formatter function will be presented in a 
late** section. 
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4. 4. 3.2 Control Logic 

The control logic consists of the master control logic 
(MCL) and the control logic (CL) circuits associated with the 
data processor. The MCL (Fig. 4.20) will perform the following 
functions : 

1. Generate a high frequency clock pulse for sub-processing 
time synchronization 

2. Generate clock pulses for frame transfer and output of 
the focal plane CCD detectors 

3. Generate clock pulses for S/H actuation in the data 
processor 

4. Generate clock pulses for input and output transfer of 
the background memory 

5. Generate clock pulses to synchronize the data buffer 

6. Generate clock pulses for data formatter 

7. Interface with the ground-based threshold detector 
adjustment input. 

8. Generate location labels to data buffer 



FIGURE 4.20 MASTER CONTROL LOGIC 
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The control logic circuits provide buffers for all the clock 
pulses and signals from the MCL to the appropriate destination. 
The MCL and control logic can be implemented using low power QMOS 
devices. 

4. 4. 3. 3 Data Formatter 

The data buffer functions as a temporary storage of the 
lightning intensity and location data. The Data Formatter (DAF) 
(Fig. 4.21) will extract the data from all the data buffers. 
Due to the low lightning occurrence rate, the data formatter has 
a relatively low data rate requirement. It will have the frame 
sampling time (milliseconds) to extract data out from the 
buffers. A time tag for each frame sampling period will be 
generated by the MCL and latched into the DAF memory. 


TIME TAG FROM | N 

MASTER CONTROL *'1 



FIGURE 4.21 DATA FORMATTER FUNCTIONAL BLOCK DIAGRAM 


64 








65 


jj 

or P'-JOh yv A'J7 Y 

into the DAF memory from the 
y this time tag. The segment 
tched into the DAF along with 
•nting 16 channels will be 

imbers. Therefore, the total 
ation identification is equal 



INTENSITY 

LOCATION 

CHANNEL k 2 
INTENSITY ATM 2 

LOCATION M 2 

INTENSITY AT 

J LOCATION Mj 

CHANNEL k« 

1 

| TIME TAG OF FRAME n 


EMORY STORAGE FORMAT 




The resolution of the time tag will be equivalent to the 
frame rate, which is 4 msec as indicated in the requirement 
section. The time tag for a 12 hour counter at 4 msec resolution 
will be equivalent to 10.8 x 10^ counts or 24 bits. The DAF is 
interfacing with the down-link telemetry. At some time interval 
the DAF will perform a data dump to the telemetry link, and the 
cycle will repeat as commanded. During a high lightning 
occurance interval, the data dump will speed up to avoid 
overflowing the DAF memory. The byte counter will perform this 
function by signaling the downlink interface when the DAF memory 
is filling up to a threshold level prior to the programmed data 
dump time. 

The total capacity of the DAF memory will depend on the data 
dump rate and the lightning occurance rate. However, since the 
occurance rate is low the memory is not required to be a high 
speed device. Consequently, it is possible to use CMOS type RAM 
to conserve power. 


4. 4. 3.4 Power Control 

The power control circuit converts the bus raw power to 
secondary power sources for the electronics. Due to the low 
power allocation it is necessary to maximize the efficiency of 
the power control circuit. Commercial DC-DC converters qualified 
for flight generally have low efficiency (50 percent to 65 
percent). Past experience with flight systems similar to LITMAP 
indicates that it is possible to design converters with more than 
85 percent efficiency and will also provide fusing and transient 
protection for proper operation of the system. 

4.4.4 Alternate Design Approaches 

There are various design approaches based on the proposed 
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system described earlier that can offer some power and size 
advantages. These approaches, however, will require further 
study and analysis when more detailed and specific system 
requirements are specified. The following paragraphs summarize 
these approaches: 

1. The effective data rate bandwidth and power requirement 
of the post-subtraction processing circuitry can be 
reduced by sharing their work load among themselves. 
All the outputs from the subtraction amplifiers can be 
fed into a multiplexing circuitry with proper sample and 
hold function and routed to several post-subtraction 
processing circuits. This approach will allow the 
reduction of the number of post-subtraction circuits as 
well as reducing their effective bandwidth requirement. 
This approach is feasible due to the low lightning 
occurrence rate. The trade-off is in the complexity of 
the multiplexing circuitry, the timing control and the 
location identification method. 

2. A variation and subset of the above method is in 
multiplexing after the logarithmic compression to reduce 
the number of A/D converters. This method will allow 
the use of an A/D converter with higher resolution. 
This approach tends to increase resolution at the 
expense of reducing data rate handling capability. 
However, at low data rate this approach seems to offer 
advantages in lower part counts and size without 
compromising capability. 

3. The background averaging technique can be treated as an 
option to provide improved system SNR and performance. 
There are alternate methods of arranging the various 
functional blocks to reduce part count and consequently 
reducing the power. 

4. Methods to reduce average power consumption of the 
system can be implemented by operating the system during 
severe thunderstorm weather, or selected seasons and 
hours of the day only. The system can then be commanded 
to a power-off or standby mode during the quiescent and 
idling period. 

4.4.5 Other Design Considerations 

There are other design considerations for the LITMAP system 
in addition to the approach and the hardware realization 
discussed above. 

4. 4. 5.1 Product Design 

The packaging of this system is Important in providing a 
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sound mechanical support and thermal as well as radiation 
protection. TRW has wide experience in designing flight systems 
operating under similar stringent requirements and as such the 
design approach will be given the most careful consideration. 

4. 4. 5. 2 Spacecraft Interface 

The LITMAP system as a payload will be interfacing with the 
various capabilities of the bus. The interface design and system 
management will be an important design and development task. 
Detailed implementation and design will be carried out when 
specific design requirements are established. 
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4.5 RADIATION EFFECTS IN SILICON DEVICES 


4.5.1 General 

The synchronous orbit radiation environment can produce a 
number of effects/defects within silicon based components which 
can broadly be grouped into two categories: 

CATEGORY PRIMARY EFFECT SOURCE RESULT 

soft electronic ionization e, p, gamma free electrons 

hard atomic collision fast neutrons displaced atoms 

References such as Chapter 6 of Barbe (1980; discuss these 
effects. Some brief comments can be made as follows. 

4.5.2 Soft Damage 

The soft damage is often transitory, but under certain 
conditions can result in a cumulative effect, e.g., trapped 
charge accumulation at a dielectric interface. Positive bias and 
surface devices are most susceptible to this trapped charge 
effect. Designs using negative bias and buried channels help 
minimize this effect. Gamma rays and a few MeV electrons can 
also sometimes cause hard damage referred to as point defects. 
The point defects usually are very limited in extent and are not 
of a major concern. 

4.5.3 Transient Effects 

A transient effect that must be considered is that of a soft 
e/ent producing enough free charge within a given pixel which is 
accumulated and interpreted as a lightning event. A radiation 
dose y at which the CCD pixel well will become saturated can be 
estimated by 

y = N FW/fty(L + W) A 

where 

Npw s full well capacity in number of electrons 
g s electron/hole pair (ehp) generation rate 

A s surface area of pixel accumulation site 
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W = depth of pixel accumulation site 

L = diffusion length 

For a buried channel device with Np W = 5 x 10^ electrons, L 

= 50 um, W = 12 um, A = 90 um 2 and g = 4.3 x 10^ ehp cm~3 rad 

(Si)” 1 , a well is predicted to saturate after a dose of 0.2 rad 
(Si). There are about 10^ seconds in a five year period. For a 
1o6 rad (Si) accumulated do3e, this gives about 10 ”2 rad/sec, 
about 4 x 10“5 rad/frame or about 6 x 10 -11 rad/pixel-frame. 
Although a single alpha particle could generate 10& ehp or a 
cosmic ray event (> 100 MeV) around 105 ehp, the event rates are 
so low as to be lost in the nominal background induced false 
alarm rate. 

The CCD recovery time Tr after a transient pulse of ionizing 
radiation can be approximated by 

T R = 8 y V coll /N Fw f c 

where V C oll * s t ^ ie we H collection volume. For a 10? Hz clock 
rate, Tr z 10"? sec. Actual measurements on devices has 
Indicated that recovery times vary from less the 0.1 milliseconds 
to over a second. Device design and selection (c.f., SIRTF, 
1978) seems quite reasonable to hold the recovery time to below a 
single frame time for these soft effects. 

4.5.4 Hard Damage 

Hard damage is long lasting and results due to fast neutron 
collisions within the bulk material causing a displacement of 
atoms from their normal positions in the semiconductor crystal. 
This displacement damage leads to significant decreases in 
carrier concentration, carrier mobility and minority carrier 
lifetime. These effects become significant for greater than 
1011, 1 0 1 3 , and 1 0 "* 5 neutrons/cm 2 respectively where about 50% of 
their energy goes into these "cluster (displacement) defects". 

4.5.5 LITMAP Electronics Hardening Approaches 

Two approaches can be taken in order to tolerate the 
radiation environment. One is to harden the devices so they can 
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withstand the exposure and the 3econd is to shield the devices so 
that the effective exposure levels are dropped to a tolerable 
range . 

Figure 4.23 (Space Systems and Technology Workshop II, 1982) 
shows the trapped particle dose as a function of aluminum 
shielding thickness with synchronous orbit solar particles 
contributing an additional few thousand rads. CMOS 
microprocessors, which are currently planned for future missions 
are able to withstand up to approximately 1 0 5 rads. Local 
shielding of 150 mils of aluminum will reduce the exposure to 
below this level, hence permitting 5-10 year missions. 

Currently available unhardened CCD's with high sensitivity 
and low noise are available which can perform up to about 10^ 
rads radiation exposure. Localized shielding of 300 mils of 
aluminum will reduce the expected dose to below this level. This 
will allow required performance for missions in the 5-10 year 
range . 

II and Hughes have been developing radiation hardened CCD 
structures. Presently Hughes has some multiphase, buried channel 
devices which perform fairly we- 1 up to 3.5 x 1 0 5 ra d (Si) 
(Chang, 1980). The TI virtual phase devices have been tested 
(McGartn, 1981) and performed fairly well up to 10^ rad (Si). If 
these hardened devices are available in LITMAP configurations, 
then the shielding requirements can be relaxed. For the 
synchronous equatorial orbit application, the spacecraft material 
plus a little shielding should provide the necessary reduction in 
radiation to permit five to ten year lifetimes. Hence, although 
the radiation environment must be considered, it does not appear 
to be a serious problem for the LITMAP application which does not 
have to survive a low, high inclination earth orbit or a 
thermonuclear event. 
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5.0 IMAGING CHARGE-COUPLED DEVICE TECHNOLOGY ASSESSMENT 


In considering suitable types of imaging charge-coupled 
devices for the Lightning Mapper application, discussions were 
held with domestic manufacturers to determine the applicability 
of currently available devices and also to determine the 
possibility of developing custom devices which would be optimized 
for this application. 

Discussions were held with vne following respresentati ves of 
the manuf acturers listed: 

a) Dr. Dean Collin3, Texan Instruments, Central Research 
Laboratory 

b) Dr. E. D. Savoy, RCA Electro Optics and Devices Solid 
State Division. 

c) Dr. Garry W. Hughes, RCA David Sarnoff Research 
Laboratories 

d) Dr. David A. Robinson, Hughes Aircraft Co., Industrial 
Electronics Group Technology Center. 

e) Dr. Rudolph Dyck, Fairchild CCD Imaging 
5.1 CCD REQUIREMENTS 

The primary requirements of the CCD focal plane array for 
the Lightning Mapper sensor application are: 

a) Large format, to obtain wide angular coverage with high 
spatial resolution. 

b) Frame transfer configuration, to enable rapid transfer 
of induced electron charge from the imaging area of the 
CCD and thus almost continuous integration of radiation 
from the observed scene. 

c) High quantum efficiency in the near-infrared range c f 
the spectrum, to permit observation of energy in either 
the 777 nm O 2 band or the 868 nm N 2 band. 

d) A high frame rate, in the order of 3 to 5 msec, to 
permit an optical time history of the lightning flashes 
to be obtained. 

e) Partitioning of the CCD to use multiple parallel output 
ports and preamplifiers to enable rapid readout of the 
on-chip storage section and high frame rate. 

f) A large pixel size, with large full well (electron 
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charge) capacity, permitting use of a large optical 
system for a given angular subtense of the pixel. The 
combination of the large optic and high full well 
capacity results in a high S/N ratio. 

5.2 SENSOR CONCEPTS 

As previously discussed in Section 4.3, the CCD focal plane 
array and sensor concepts which have been developed in this study 
consist of the following (Figure 5-1): 

a) Two Large Format CCD ' S 

Two large format (800 X 800) pixel imaging CCD’s 
modified to p r ovide on-chip frame storage and multiple 
output ports and preamplifiers to enable fast 
(parallel) readout and thus a very high frame rate. 

b) Four CCD ' s with One or Four Optic:’. 

Four individual CCD’s modified by partitioning to 
obtain multiple output ports and high frame rate. The 
optical fields of view might be combined by optical 
beam splitting, or individual optical systems could be 
used for each of the four CCD's. The former approach 
is difficult to accomplish using an optic with a low f « 
number and presents problems in design of the narrow- 
band spectral filters due to wide field of view of the 
4-chip array. The latter approach eliminates optical 
beam splitting and results in a reasonable angular 
fields j.‘ '>iew for the narrow-band spectral filters. 

c) Four - «ith I ntermittent Integration Sampling 

Four RCA frame transfer imaging CCD's, in either of two 
sensor configurations as discussed in b) (above), in 
which intermittent sampling of the incident optical 
radiation would be accomplished by the use of a short 
integration interval in conjunction with a longer 
readout time of the on-chip storage section. The short 
exposure time would be obtained by dumping most of the 
induced charge in the image section into a guard ring 
structure above the image section of the CCD, and using 
the remainder of the normal exposure interval for 
integration in the conventional manner. Thus, with 
operation at a frame rate of 60 Hz, a useable 
integration time of 4 msec out of the normal 17 msec 
Integration tire could be obtained. 

d) One or Two Small CCD's With Small FOV 

One or two CCD's with a single optical system covering 
only a smell field of view. Here currently available 
CCD's would be used, but only very limited areal 
coverage would be obtained with a mechanical storm 
tracking mode being appropriate. 
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NARROW 

BAND 

FILTER 


ORIGINAL PAGE IS 



CCD #1 

4, 8, OR 16 OUTPUT TORTS 

400 x 800 PIXEL STORAGE 
SECTION 

400 x 800 PIXEL IMAGE 
SECTION 

CCD #2 

400 x 800 PIXEL IMAGE 
SECTION 

400 x 800 PIXEL STORAGE 
SECTION 

4, 8, OR 16 OUTPUT PORTS 


A) TWO LARGE FORMAT CCD's 



B) FOUR CCD's WITH EITHER A COMMON OPTICAL SYSTEM AND SPECTRAL 
FILTER OR INDIVIDUAL OPTICS AND SPECTRAL FILTERS 

C) SAME AS B) BUT WITH SHORT EXPOSURE TIME OF CCD 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



D) ONE SMALL CCD WITH A NARROW FIELD OF VIEW 
GIMBALLtD ON TWO AXES 


FIGURE 5.1 SENSOR CONCEPTS 


5.3 CURRENTLY AVAILABLE CHARGE-COUPLED DEVICES 

Table 5.1 lists the characteristics of currently available 
imaging CCD's which can be considered for the Lightning Mapper 
application. These are indicative of the technological 
capability of the manufacturers listed in terms of CCD 
configuration and performance. Some of the devices can be 
directly applied to the LITMAP application, while others would 
require modification. Each device in Table 5.1 will be briefly 
discussed. 

The largest format in a single CCD is that of the 800 X 800 
pixel device from Texas Instruments. The device is in a line 
transfer configuration, and modification to a frame transfer 
configuration with multiple readout ports : ^uld be required. TI 
is also developing a i 024 X 1024 device for NASA, but this is 
still in the early developmental stage. A disadvantage of this 
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TABLE 5.1 Summary - Characteristics of Available Imaging CCD 


ORIGINAL PAGE JS 
OF POOR QUALITY 





■ MOST SUITEO FOR lITMAP APPLICATIONS. MODIFICATION REQUIRED 

A DEFICIENT FOR LITHAP APPLICATION (SEE APPENDIX A. B AND SECTION 2.3) 1. Q.E ► SENSITIVITY 

A OPTICAL THINNING/COATING, ELECTRODE CHANGE (SnO~). CYLINDRICAL LENSE CONCENTRATION 2 FW ► OYNAMIC RANGE 

COULD BE LARGER 2 3. PIXEL SIZE ► APERTURE 

O UPGRADE WITH NEW OUTPUT AMP MITH SMALLER CAPACITANCE FOR HIGHER SPEED/LOWER NOISE READOUT 4 FREQUENCY READOUT TIME 



CCD is the small pixel size, limited by the size of a device that 
can be produced with reasonable yield. Yield is limited by the 
degree of purity (number of defects) in the silicon base 
material . 

Two CCD's are in development by TI for use in commerical 
television, the TC 201 (for U.S. standards) and the TC 202 (for 
European standards). The primary difference is that the U.S. 
version does not use blooming control, eliminating the need for 
anti-bloom drains within the imaging format, resulting in a more 
contiguous imaging area (higher packing density) and increasing 
the effective quantum efficiency. Pixel sizes of the two devices 
are comparable. 

The SID 501D has recently been announced by RCA, replacing 
the previous SID 53612. The former has an overlapping gate 
structure and improved preamplifier, while the latter has a 
planar gate structure. This device is listed primarily to 
indicate the current capability of RCA in this field. This 
device uses a thinned substrate (laminated to glass and backside 
illumination to obtain high quantum efficiency. This sometimes 
causes problems in the near-infrared range of the spectrum due to 
interference fringe patterns occurring in both the glass laminate 
and silicon which result in high fixed pattern noise. Anti- 
reflection coatings and adjustment of substrate thickness may 
help resolve these problems. For the LITMAP application, 
frontside illumination and thick CCD substrate maybe preferable 
to eliminate the problems associated with deeply penetrating 
infrared radiation. RCA has previously manufactured such devices 
for use in CCTV cameras. 

The salient characteristic of the Hughes H4068 CCD is that 
of radiation hardness, being processed to withstand both natural 
and weapon-induced radiation. This device has been developed for 
use by TRW in star trackers used for satellite attitude control 
under the MADAN contract to the USAF Space Division. The 
limitation of this device is that of the line transfer 
configuration, without on-chip data storage. However, it is 
feasible to add on-chip frame storage to this type of device. 
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The Westinghouse type 4068 was developed under NASA funding 
in 1977. Although small in format, the device has the advantage 
of high quantum efficiency due to the use of SnO gates, which are 
very transparent in the visual range of the spectrum. This 
device would be applicable for the LITMAP sensor configuration 
where only a limited optical field of view is used. 

Fairchild is well advanced in development of the CCD 222. 
This device is used in both commercial CCTV cameras as well as m 
USAF military aircraft and helicopters. Considerable Navy 
funding has resulted in the development of integrated CCD logic 
and driver modules, as well as modules for video processing, 
resulting in versatility in available camera configurations. The 
primary limitation of the CCD for the LITMAP application is the 
use of opaque interline frame storage registers in the image 
section, which block out 50% of the imaging area. However, this 
limitation might be overcome by the use of a faceplate consisting 
of a series of miniature cylindrical lenses precisely registered 
to the CCD imaging area, similar to the method suggested by Cross 
et al (1982) 

5.4 DISCUSSIONS WITH MANUFACTURERS REGARDING DEVICE OPTIMIZATION 
5.4.1 Texas Instruments 

The present 800 X 800 Galileo Project CCD uses an amplifier 
which is optimized for a slow readout rate (50 KHz) with the CCD 
at low temperature to reduce dark current. For the LITMAP 
application the preamplifier would be replaced with a wideband 
preamplifier of the type used on the TC 201 and TC 202 devices. 

Partitioning the device for parallel readout with multiple 
preamplifiers presents no problem. Four, 8 or as many as 16 
output ports could be used. Clocking can be performed at rates 
as high as 17 MHz without loss of charge transfer efficiency. 
However the use of correlated double sampling at this high 
readout rate is not feasible, as both sampling and transfer 
pulses must occur within the pixel readout interval. Without 
double sampling, a noise level in the order of 125 electrons 
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rms/sample should be realized with a readout of 6 MHz. 

TI would want to maintain the small (0.6) mil pixel size, as 
the large format results in a large chiD ( 2 - ^ cm 2 including I/O 
provisions) and yield is limited by chip ize. Some reduction in 
sprtial resolution would be expected due to the roll-off of the 
modulation transfer function in the near infrared range of the 
spectrum due to small pixel size. This would be more critical at 
the 868 nm wavelength than at the 777 nm wavelength. 

For a CCD of this large size the cost of a new set of 
photomasks is in the order to $100K, and t^tal program cost would 
be in $1M category. The duration of a development program of 
this type would be 1 to 2 years, depending upon requirements. 

Interest of TI management in committing resources to the 
program would be based upon validity as a business proposition. 
TI would appreciate sharing in development of the total sensor 
design concept, rather than just accepting procurement 
specifications for CCD development, as previous related 
experience might be of benefit to the program. 

Regarding t. e TC 201 and TC 202 types of devices, emphasis 
is being placed upon initiating production for use in commercial 
CCTV applications. Devices are available at relatively low cost 
with short delivery times. For this application the TC 201 would 
be preferred due to higher packing factor and quantum efficiency. 
These devices could be used in the four-chip focal plane array 
configuration with operation in a frame transfer mode if somewhat 
lower resolution is acceptable. Total format for four CCD’s 
would be 648 X 490 pixels. The TC 202 would provide a total 
format of 780 X 584 for four CCD’s 

5.4.2. RCA Electro Optics and Devices Solid State Division 

RCA has discontinued work on all frontside illuminated CCD's 
at the Lancaster, Pa. facility. All devices are thinned and 
backside illuminated. Major products are the SID 504 for CCTV 
applications, with small pixels (20 x 20 urn) and the SID 501 with 
large pixels (30 x 30 urn) for scientific applications. Both 
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would require special attention for this application due to the 
use of backside illumination with a very thin (10 um) substrate 
which results in severe fringing when observing near-infrared 
radiation. However, Dr. Savoy's group is considering development 
of a thicker scientific CCD wi' h a substrate thickness of 30 um. 
In combination with a higher resistivity material, this type of 
device might be more suitable for use with red or near-IR 
radiation . 

For this operation, Dr. Savoy suggested the following type 
of operation for the LITMAP application. The image and storage 
sections and output register of RCA CCD's are surrounded by a 
guard ring shaped like a picture frame which can be used as a 
drain to absorb unwanted photoelectrons from the image section. 
Assuming a normal integration time of 17 msec with a 60 fps field 
rate, the image section could be allowed to integrate optical 
radiation for 13 msec. Presumably a large amount of charge would 
be accummulated due to solar illumination incident on the earth. 
This could then be rapidly clocked upward into the guard ring 
(drain) and absorbed. The remaining 4 msec could be processed to 
obtain measurements of storm activity, and the induced electron 
charge in the image section would be clocked into the on-chip 
storage section and read out slowly in the usual frame processing 
manner. This method of operation would permit use of an on-the- 
shelf device in early feasibility tests of he LITMAP sensor. 

In regard to accepting custom work in CCD development, Dr. 
Savoy stated that management at RCA is placing extreme emphasis 
on getting the SID 501 and SID 504 devices into production, and 
that the Lancaster facility could not accept custom work for the 
foreseeable future. This situation would probably change in a 
year or two. 

5.4.3 RCA David Sarnoff Research Laboratories 

Dr. Gary W. Hughes was quite familiar with the requirements 
for the LITMAP sensor, having previously worked with Dr. Wolfe of 
the University of A. izona on this application. 

He agreed with the concept of using frame transfer CCD's 
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with multiple output ports. The parallel transfer rate (from 
image to storage sections) is limited to a maximum of 800 KHz by 
CCD capacitance and driver capacity. Output noise would be in 
the order of 200 electrons rms without correlated double sampling 
and 35 to 40 electrons rms with double sampling at a readout rate 
of 7.5 MHz. It may be possible to use double sampling at readout 
rates as high as 10 MHz. 

He has been involved in the design and development of the 
SID 501 and 504 devices, currently in initial production at RCA 
Lancaster. These devices are representative of the largest 
format which has been developed by RCA. 

If RCA were involved in development of a device for the 
LITMAP application, design would be performed at the ••esearch 
laboratories and fabrication would be accomplished at the 
Lancaster facility. 

A rough order of magnitude estimate of cost of a development 
program of this type would be $700K to $1M, with schedule in the 
order of one year. Interest of RCA in bidding on a CCD 
development program would be based upon business potential. 

5.4.4 Hughes Aircraft Co. , Industrial Electronics Technology Center 

Hughes Aircraft Co. has been involved in the development of 
radiation hard CCD's since 1975 under contract to NRL, DNA, and 
the Special Projects Office of the Navy. To date imaging CCD's 
have been developed for both the Navy and TRW (under the USAF 
MADAN project) with formats of 324 x 324 pixels in a line 
transfer configuration. 

Hughes proposes development of a device with a 400 x 800 
pixel format, with half of the area being used for on-chip frame 
storage. The storage section would have eight output ports, each 
operating at a rate of 5 MHz, resulting in a frame time of 4 
msec. At this readout rate the noise level would be no larger 
than 500 electrons rms. As with the MADAN devices, the full well 
capacity would be at least 10$ electrons. Pixel size would be 20 
x 20 micrometers and quantum efficiency at 868 nm would b<~ 30 
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percent or higher. 

A budgetary and planning cost estimate was submitted by 
Hughes to TRW for development of both imaging CCD's and 
breadboard electronics. 

5.4.5 Westinghouse Advanced Technology Laboratory 

Westinghouse is currently under subcontract to TRW for 
development of an improved version of the 5040 imaging CCD. 
Details of the device are presented in Table 5.1. This device 
would be suitable for use in the sensor concepts described in 
Section 5.2 d. In conversation with Dr. James A. Hall, it was 
determined that Westinghouse has not developed CCD’s with very 
large formats, although they would be interested in reviewing 
procurement specifications for such a program. 

5.4.6 Fairchild CCD Imaging 

The Fairchild 222 has a favorable frame transfer time due to 
its interline transfer and cylindrical lens arrays are being used 
to increase their detection quantum efficiency. The devices are 
used for a number of DoD applications. This device would be 
suitable for use in the sensor concepts in Section 5.2c but 
discussions for the development of a special device for LITMAP 
with multiple ports have not been made. The device is being 
considered for inhojse high frequency evaluation of the LITMAP 
frame substraction concept. 

5.5 CONCLUSIONS 

Reviewing the requirements defined for the CCD in Section 
5.1, we reach the following conclusions: 

a) For m at Siz e - The largest format which has been 
developed to date is the 800 x 800 pixel device from 
Texas Instruments. However, in order to maintain a 
reasonable size of the chip to obtain an acceptable 
yield in wafer processing, the pixel size is small (0.6 
x 0.6 mils). This results in the following: 

o Small pixel size reduces the full well capacity, 

which due to the use of the virtual phase 
structure is 109»000 electrons for this device 
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For a given angular field of view for each 
pixel, a smaller pixel results in a shorter 
optical focal length and a smaller lens 
diameter. As less optical radiation is 
collected, a reduction in S/N ratio results 

o The smaller pixel size results in higher optical 

cross-talk when observing radiation in the near- 
infrared range of the spectrum. Due to the low 
absorption coefficient of the silicon, the near- 
IR photons generate carriers below the depletion 
region of the CCD, which migrate laterally to 
adjacent pixels. This results in loss of 
resolution which will be evident as broadening 
of point source images. 

Thus the combination of small pixel size, low full well 
capacity, and increased optical cross-talk do not make a 
large format virtual phase CCD attractive for this 
application. None of the manufacturers surveyed were 
willing to consider development of a CCD with both a 
very large format and very large pixel size. 

On the other hand, the TI TC 201 CCD looks attractive for 
this application due to large pixel size, large full well 
capacity, and high quantum efficiency. Format size when used in 
the frame storage mode is 328 x 245 pixels. This device will be 
applicable in a sensor configuration where multiple chips are 
used . 

b) Frame Transfer Configuration, Partitioning, ancl High 
Fra me Rate - All of the manufacturers agreed that in 
order to operate in a frame transfer mode at a high 
frame rate, partitioning of the on-chip storage section 
of the CCD would be required. In the case of the 800 x 
800 pixel TI CCD, partitioning into 4, 6, " * as many as 
16 sections with a corresponding number of output 
preamplifiers was suggested. 

c) Recommended Focal Plane Array Configuration - Only Texas 
Instruments has demonstrated the capability of producing 
very large CCD's (800 x 800 pixels). However, tnese 
have disadvantages as discussed in item (a) above. TI, 
RCA, and Fairchild have demonstrated the capability of 
producing devices with a format in the order of 250 x 
300 to 400 x 500 pixels, with pixel sizes ranging from 
12 x 18 to 30 x 3u micrometers. A practical approach to 
the focal plane array for the LITMAP sensor would be to 
use a mosaic of chips of this type with either a common 
or independent optical system. 
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A. PARAMETER SPECIFICATION 


A.1 INTRODUCTION 

This appendix contains some details and rational for 
selecting the parameters for the baseline concept. The system 
analysis as summarized in Section 3 and Appendix B involved 
performing parametric trades about these baseline specifications. 

A summary of the results of this appendix is given in Section 2. 

A. 2 FIELD OF VIEW: FOV 

One of the first parameters to be established is the field 
of view (FOV). Full disc (FD) and continental United States 
(CONUS) are given in Table 2.1 as design goals with emphasis 
being placed on FD. Geosychronous orbit has also been specified 
in the SOW. From this vantage point (approximately 5.56 earth 
radii altitude) the absolute viewed disc will be less that the 
full earth diameter. This geometric viewing coverage is further 
restricted by the increasingly oblique line of sight (LOS) angle 
at the outer edges of the viewed region which is due primarily to 
the earths curvature. These factors will now be quantified. 

The earth latitude (6), associated sensor viewing angle ( <P 
) and distances in a projection plane (tan = tangent at the 
equator) are summarized in Figure A.1 with representative 
numerical values given in Table A.1. 

For equal sensor pixel size over the full <'<-■. the 
incremental observation angleAVi ,is constant. As <p iv ea> • 
the surface distance (S) (hence area' increases rapidly 
illustrated in Figure A.2A. This normalized incr* , r -■ j 
illustrated in Figure A. 3. For const .it longitude, orig. f^ng 
at the sub-satellite point, the area observed by a pixel will 
vary approximately linearly with the incremented change in the 
distance along the earth's surface, AS . Along a diagonal, both 
longitude and latitude are changing and the area increases in 
worse case as approximately AS?, This latter change is 
represented by the DA/DA 0 curve in Figure A.3« The relationship 
between ? and e is shown in Figure A.4 which also 
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FIGURE A. 2 Distance Along E 
Angle (A) and Ea 
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3 Linear and Area Variations (for equal steps) Normalized to the 
Respective Dimension at the Subsatelllte Point Relative to Sensor 
Observation Angle (A) and Earth's Latitude (B). 
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FIGURE A. 4 Variation of Earth's Latitude (fixed longitude) with Observation 
Angle. 

illustrates the rapid increase in earth's latitude as the 
observation angle increases. 

The selection of the full FOV limits must be based on 
several factors which include such consideration as: desired 

surface coverage, SNR degradation and circular earth cross- 
section vs square sensor array. As seen in Figure A.1A, if the 
full disc is covered , a latitude coverage to approximately 81 
degrees can be provided. For a nominal square sensor array, this 
choice would result in a significant number of pixels which will 
extend beyond the earths surface, hence not contributing to the 
primary objective of lightning detection. Although circular 
arrays may not be out of the question, they would require 
complete new design with special processing circuits. Since, as 
seen in Figure A.3B, the high latitude areas will have rapid SNR 
degration due to the observed area Increase coupled with finite 
source size, a smaller FOV coverage would be a more practical 
choice. 
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If a square inscribed within the full disc is chosen, the 
coverage would be limited to about 44 N (the east-west can always 
be covered by an increase in number of satellites). This 
excludes a portion of the CONUS. Figure A.3B shows that the 
constant longitude linear dimension (hence area in this case) 
increases by a factor of 2 at about 52 degrees . The upper limit 
of the CONUS is at about 49°N. The selection of 50°N would by a 
reasonable compromise amoung the various factors mentioned 
previously. 

FOV Selected : (half angle) — 7.37°(50°N latitude) 

A. 3 AREA PER PIXEL: A B 

From Table 2.1 a nominal ground resolution goal of 5km is 
specified which for a 5C°N FOV would result in a linear pixel 
count, N, of about N=S/AS N =5770km/ 5km = 1,114. The full disc 
FOV would require an array with twice this dimension of 2228X 
2228. As shown in the technology assessment section, single, 
double or quad arrays satisying this criteria do not exist at 
this time. 

Special sensor developments for DoD applications are 
directed at the 1250 square array catagory, some of which may be 
amenable to LITMAP requirements (submicron wavelengths). These 
(e.g., RCA PtSi and Pd2 Si technology) show promise (due to 
restriced image spreading) for furture generations of LITMAP 
(post 1990). 

Within the present technology, the most practical single 
array size would be a 400 X 400 (maybe 500 X 400) active sensor 
array with matching on chip frame storage. Chips with 390 X 292 
active region with frame storage and moderate pixel size do exist 
which would require only the addition of parallel output ports to 
handle the LITMAP framing time. 

Based on the technology assessment (see Section 5), a square 
arrangement of four 400 X 400 units was selected for the system 
trades. The average resolution for this array size is then given 
by the distance/number s S/N = 5570/400 s 13.9 km/pixel. Using 
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Figure 2.4 this nominal pixel size translates into a as 
( subsatellite resolution) of 12.2 km and 10.3 km when the nominal 
13-9 km resolution element is located at the center of the upper 
hemispheric FOV (25°N) or at the center of CONUS (37.5°N) 
respectively. 


N- ARRAY SIZE SELECTED: 

N = 400 

(North) plus 

frame storage 

Sub- 

Satellite 

Central 

CONUS 

50°due north 
Sub Satellite 
Point 

Scaling 

1.00 

1.35 

1.32 

Resolution Selected (km) 

1 1 

14.9 

2C 

Area/p: xel (km^) 

121 

163 

220 

A. 4 SOURCE CHARACTERISTICS 




The key source properties in the SNR relationship are: 
source area, lightning pulse length, energy distribution, 
detection wavelength, associated energy per photon and 
source/pixel fill factor. Each of these will now be quantified. 

A . 4 . 1 Source Area: A s 

The source area specified in Table 2.1 is that of the source 
illumination filling a representative cloud cell, given as 5km 
radius. The ideal situation is that the source fill the full 
pixel, that is keep the background area Ajj as small as possible 
without dividing an lighting event between f'o or more pixels. 
This consideration is important when defining the true 
statistical threshold detection capability. For larger lightning 
events signifcantly above threshold, subdivided events can be 
post-detection summed to give the original even energy. 

Data on effective source area is very limited. One approach 
is to assume that the source always fills the pixel FOV, but as 
the FOV increases, this assumption could cause significant error 
in assessing a design capability. Lopex (1977) has given some 
data as shown in Figure A. 5. Fifty precent points* range from 
below 0.1 km to 25 km for various parts of the world with most 
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falling in the 2 to 10 km range. For baseline systems analysis, 
the filled pixel condition was assumed and a parametric analysis 
made about this baseline. Due to pixel geometry and considering 
the source size, a square geometry approximation was used. 

A s = nominal source area: 121 km^ at sub-satellite point 


A . 4 . 2 Geometric Splitting F actor : v 

In the ideal situation, the event is assumed to be centered 
on the detection pixel. In the lightning mapper application, 
the centroid of an event can fall randonly, anywhere with the 
sensing pixel area. Even if the event's area is smaller than the 
pixel's footprint at the source, there is still the probability 
that the event will straddle the boundary between two or more 
pixels. Streuber and Bailis (1980) has addressed this situation 
with the quantitative results shown in Figure A.6A. When the 
pixel size g is equal to the source size 6 S ( <5p/<5s = 1 where 


Jult J977 


RAUL ERLANDO LOPEZ 




Accumulated frequency distribution) of cloud and echo 
horizontal dimension. The straight lines correspond to the log 
normal distributions that beat it the different data sets. 


FIGURE A. 5 Distribution of Horizontal Cloud Dimensions. 
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the source size has been choosen as the four sigma point) there 
will be on the average, about 62 t of the source energy falling 
within a specific pixel bounds. The implication of Figure A.6A 
is that the larger the background area, the better the results, a 
condition which would be somewhat true for night time conditions. 

But under daytime, solar background conditions, a larger 
background FOV means more noise, hence a decrease in the SNR. 
This overall factor is qualitatively illustrated in Figure A.6B 
for the situation where solar background noise predominates over 
all other noise sources. For background conditions where 
>> 1 , the SNR does decrease. As the background FOV area is 

decreased however, there is a relative maximum in the SNR/63 
factor, after which the SNR decreased slightly until SNR/tfg 
becomes, in the limit, approximately 1/2 * . This means, for a 
fixed source size and as the background area becomes 
significantly smaller than the source area, there is a point 
where the pixel will always be filled with equal area of 
background and the lightning source at which time the SNR will 
converge to a fixed value. In practice, noise sources other than 
background would result in a lower limit for the resolution size 
below which no SNR improvement will occur. In addition, sensor 
technology constrains the selection of the array pixel count. 

As disussed previously the source dimension is likely to be 
smaller than the background FOV for attainable array 
configurations. For sub-satellite fill factors less than one, 
the 83/63 3ca l es proportionally (for constant longitude). This 
fill factor correction is an important parameter in the 
parametric trades. 



6 Blh 

V 

LOCATION 


0.63 

0.«M 

optimum 

Selected 

1 

0.62 

at sub-satillite point 

Values: 

1.35 

0.70 

at 37.5 N 


1.82 

0.78 

at 50 N 
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A. 4. 3 Lightning Pulse Duration; !^ 

Lightning pulse duration contributes to the definition of 
the shortest integration time that can be used if time splitting 
of the pulse is to be avoided or at best optimized. Figure A. 7 
summarizes the effective pulse width of lightning for some recent 
U-2 data (Gillapsy, 1/83) and for data from Turman (1976). Both 
imply an increase in pulse duration with an increase in total 
energy. Care must be excersized here, because the peak energy 
and peak power data were not necessarily from the same event. 
Subsequent data from NASA/ MSFC (Gillaspy 6/83) reproduced in 
Figure A.8A for maximum stroke per flash and Figure A.8B for all 
strokes indicates an interesting alternative. Even though there 
is a variance of effective pulse duration for any energy level, a 
constant average, effective pulse duration could be specified 



FIGURE A. 7 Effective Pulse Duration (peak energy/ Integrated power) 


98 









which would apply to all maximum strokes independent of total 
energy envolved. Figure A.8B indicates that this conclusion 
could be loosly applied to all strokes. An assessment of the 
NASA provided data for the full width at half maximum indicates a 
value of 0.54 millisconds at the 50U point which is relatively 
close to previous data from Turman of 0.45 milliseconds. 

T l - selected value: 0.54 milliseconds 

A. 4. 4 Pulse Splitting : £ 

Nagler (1981) addressed the fact that lighting events occure 
randomly in time relative to the frame sampling time, hence 
resulting in time splitting of strokes between two frames. The 
nearer the integration time is to the pulse duration, the more 
prorainate thjs problem becomes. As mentioned for the spatial 
splitting, the occurance of this phenomenon in most important 
near threshold where the loss of part of a pulse could result in 
a loss of an event with the concequence being an effectively 
higher threshold. Sensor dead time (e.g. during frame transfer) 
will also contribute to the loss of data. An assessment of some 
of the pulse waveforms provided by NASA/MSFC indicates a quisi- 
triangular waveform. 

The time splitting factor for a triangular wave form is 
given by: 

C • T, / T l - (T, / T l ) 2 /3 

{ . 1 - ( \ / T, ) /3 

In the SNR relationship, the full time dependency is given by or 

tnfi 7 in normalized terms C(T^/T^ )*/(T^)* . The factor Tl is 
assumed constant and the normalized { (Tl/T*) 1 ^ i a plotted in 
Figure A. 9. For maximum SNR, T^sTl* but small variances in 1 l 
would result in significant variation in pulse detectivity. The 
SNR falls off extremely rapidly for Ti<Tl* To minimize a large 
variation in SNR as the event pulse width changes, a value of 


0 < T. < 

• 1 


;oo 




FIGURE A. 9 Time Factor Effect on SM° 


Ti>>T L should be chosen, a factor greater than 5 would be 
recommended . 

Other wave forms result in some variation of the constant 
1/3 in £> but the peak will still occur in the same place(at 
Change in steepness of the curve and a change in the 
absolute magnitude of this time splitting factor would occure, 
however, with no significant change in the conclusion. 

4 - Stroke time splitting factor: (1 - Tl/ 3T i) 


A. 4. 5 Energy Strength : E a 

Table 2.1 indicates two candidate wavelengths, 3683A and 
7774A. A threshold of 107 watts (total optical power) is also 
given. NASA U-2 data (Gillaspy 1/83) indicates the following 
values near the 8683A line: 


uJ/m2-sr 

4 

13 

24 

52 

100 


Probability Greater Than 
901 
50% 
mean 
90% 

Maximum Observed 


The total optical energy values were derived by assuming that a 
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FIGURE A. 10 Energy Distribution for Lightning. 


nominal value of 15% of the optical energy is contained within 
the 8683A line region. Figure A. 10 compares ground based 
lighting data with that collected by NASA/U-2 and that from space 
analyzed by Turman. Both above cloud collected data show a 
narrower dynamic range, a fact which may be real or due to a 
limited data base. Some of the U-2 d^ta showed only a 12:1 
dynamic range, but due to the limited data set, a minimum dynamic 
range of 100:1 should be used for the system analysis. A 
threshold goal of 10? watts is given in Table 2.1. This 
translates into (15 percent of energy in the line assumed) 
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T 


L 


Milliseconds 


1 


0,5 


Source Area 

km^ 


Energy in Single Line 
w J/m^-sr 


25 

100 

25 

100 


9.5 

2.4 

4.7 

1.2 


where a flat plate radiating from both surfaces has been taken as 
a typical cloud deck with an internal lightning stroke. The 0.5 
millisecond, 5km square source area give a integrated source 
strength comparable to the U-2 measurements. This result must be 
reinforced by addition field data. Modelling by GUO and Krider 
(1982) indicates that the source may not be Lambertian, but could 
have preferential radiance in a cone about the vertical, a result 
due to multiple internal scattering. The 1982 NASA/U-2 data wi^l 
be used as the baseline value for the threshold integrated energy. 

nE s ~ Radiant Energy: 4 uJ/m^-sr NASA/U2 Data 

X - Wavelength : 8683 A 

Ephot" he/ A : 2.29 E-19 Energy/Photon 

J/photon 


A. 5 SIGNAL NOISE RATION: SNR 

The average false alarm rate FAR is defined as the average 
number of times per second the output noise electrons exceeds the 
threshold setting of the detectors. This can also be '••xpressed 
in terms of the SNR as: 


FAR 


e -SNR T 2 /2 


This FAR rate decreases rapidly as the threshold value (SNRj) is 
raised. The probability of detection Pj is approximately the 
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probability f hat the signal plus noise exceeds the threshold at 
the instand of signal peak, giving 

P d » (1/2) (1 + erf ( ( SNR s - SNR T ) / VT ) ) 

where SNRj is the threshold SNR setting and SNR$ is the SNR of a 
given signal (note that when SNR$ = SNR 7 , P<j = 0-5, i.e, 50%). 
Figure A. 1 1 graphically illustrates this relationship between FAR 
and SNR, For FAR <<1, the SNR requi. ement is not very sensitive 
to the absolute value of the pulse duration between 0.54 and 1.0 
milliseconds . 

A FAR design goal of less than 10$ is given in Table 2.1. 
A. 2 summarizes the SNR required on a pixel by pixel basis to 
provide this FAR. A nominal value of the threshold SNR is around 
6.0. Note that the smaller the area which a pixel views, the more 
stringent the SNR value due to the lower lightning activity in 



FIGURE A. 11 FAR vs. SNR. 
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EVERT RATE 


3/R 



FLASH/ka^sec 

FLASH/SEC PIXEL 1 

FOR 101 FALSE EVERTS 
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(10 ka) (20 ka) 

(10 ka) 

(5 ka) 

GLOBAL 

2 10~ 7 

2 ID’ 5 0 
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6.2 

6.0 

US 

2.3 10- 7 

2.3 10-5 

lO"* 1 

6.2 

5.9 

STORM COMPLEX 

A 10~ 5 

A «0~ 3 1.6 

10“ 2 

5.3 

5.0 

SEVERE STORM 

1.7 10“ 3 

1.7 10- 1 6.8 

10- 1 

A. 6 

A. 2 

RASA/M3FC (Hugh Christian) 






GLOBAL 

2.9 10“ 7 

2.9 1 0“5 1.2 

10“* 

6.2 

5.9 

M50/*«c total) 






that reduced area. 

The respective values 

of SNR 

for footpri 

of 5,7,10,13, and 20km are 6.4 

, 6 . 3 , 6 . 2 , 6.1 

and 6, 

0 f or 

a nomi 

global rate of 2 X 

10"7 flashes/km2 sec. 

The higher 

the lo 


activity, the lower the SNR can be set for the same FAR. 


There is also a 90% probability of detection specification. 
Figure A. 12 relates the probability of detection relative to the 
SNR threshold setting (i.e, SNR at P d = 50 %). From Figure A. 12, 
an event with a SNR greater than 45% above the thre. hold setting 
will have a -probability of detection greater than 90 a 

SNR - Baseling value: 6.0 


A. 6 BACKGROUND: Ig, ALBEDO 

The solar spectral irradiance can be found from a number of 
handbooks. Approximate values for Iq tt (W/m2 - um) 

X m=0 m=2 (zenith = 60°) 

8683 950 800 W/m2- U m 

7776 1160 907 W/m2- U m 

Most of the cloud tops are above a significant part of the 
atmospheric air mass, in particular, above much of the earths 
water vapor (c.f. Figure A. 1 3 ) - Hence a worse case but realistic 
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sole • condition would be m = 0 giving for 8683, 
l B = 302 W/m 2 -um-sr. 

The albedo (solar reflectance) varies for different 
backgrounds. Of specific interest are clouds, but a pixel may 
view a broken cloud deck. Some albedo values are 


OBJFCT 


X =7776 

8683A 

irfin . er snow/ ice 

0.85 

0.8 worse care background 

C 1 o .. d s 


0.7 

0.67 normal cloud cover 

Summer ice 


0.5 

0.45 i 

Vegetation 


0.4 

0.5 /surface background 

Soil /Rocks 


0.3 

0.33) 

A nominal 

baseline 

value of 

0.7 was used for the sys 

parametric 

analysis. 

8683A 

7774A Units 


I ( - Selected Values: 300 396 W/m 2 - um-sr 

A bedo - Selected Values: 0.7 0.7 


A. 7 OPTICS: 1 /NO. ,K 

The two primary characteristics of the optics in the SNR 
relationship are the F/No (relating focal length and aperture 
diameter) and the optical throughput efficiency, K. 

The preliminary analysis indicated that the F/No would need 
to be pushed to optimize the SNR. This is one of the few factors 
which cont: ibutes propor tionaly to the SNR rather than by the 
square root. Pushing the F/No. often means mor glass, higher 
fabrication precision and mechanical limits. A F/1 is often 
considered the lower limit with F/1. 5 being more practical, 
separate design activities have indicated that an F/0.78 system 
might be doable for the LITMAP requirements. Refractive and 
reflective alternatives are possible. The requirements of the 
filter must be considered when defining the optics. 
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The optics system will probably have many surfaces, hence 
appropriate coatings to maximize throughput must be used. 
Coatings with K=0.9S to 0.99 transmittance per surface should be 
attainable . 

F/No-Baseline : 1.0 

K-Sele.ted value: 0.98 

A. 8 FILTER: 

When considering the filter requirements, care must also be 
taken that the throughput of the filter is also considered. 
Although the optics and filter were not part of this study, their 
significant impact on the overall system trades necessitated a 
look at their impact during a related IR&D activity. 

Table A. 3 compares the overall effects of two filter types 
and several optical band passes. The interference filter must 
have a bandpass sufficient to handle the angular effects (c.f. 


TABLE A. 3 

Filter/Optics Figure 

of Merit 


riLTEA TYPE 

* 

K 


Fliure of Aerlt 

/ran 


(A) 

FILTER 

FILTER ♦ OPTICS 

INTERFERENCE 

25 

0.90 

0.57 - 0.72 

0.15 - 0.17 


5 (rectangular) 

0.70 

0.56 

(3) 0.83 

(5) 0.31 - 0.33 


2.5 (gausalan) 

0,80 

0,25 - 0.32 

0.31 - 0.35 

BIREFRINGENCE 

(SOLC) 

2 

0.25 

0.16 - 0.20 

0.28 - 0.32 


1 

0.25 

0.16 - 0.20 

0.8 . 0.88 


* of optical lint being isolated if assuaad to ba on# or two engstroas. 

*• Assuaad optics with 11 eleaents (22 surfaces); (0.99) 22 * 0.9; (0.98) 22 * 0.68 


Figure A. 14 while the SOLC filter can handle the larger FOV while 
maintaining a narrow bandpass. 
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^(DEGREES) 

FIGURE A. 14 Angular Shift for Interference Filters. 



For wider FOV's (e.g., I8°full angle) the SOLC has a 2 to 3 
advantange over the required 25A bandpass interference filter. 
For the narrower FOV (as exists when the full disc is split four 
wayswith separate optics for each, see Section 4) both filters 
provide comparable performance with the development risk being 
larger with the SOLC filter. If more than 3 A bandpass is 
required, interference filters are the most cost effective 
approach . 

For band passes of a few angstroms, a single line of the 
lightning spectra can be isolated. An approximate value of 1 5t 
of the total optical energy is associated with each of the 
candidate lines of 3683 and 7774 angstroms. 

A A - Selected : 5 A 
i? - Selected: 0. 15 
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A. 9 SENSOR CHARACTERISTICS 


A. 9.1 Quantum Efficiency; Q 

The quantum efficiency term includes the overall conversion 
of the total incoming photon flux to charge carries within the 
detector. This not only includes the absorption properties of 
the active region but also accounts for the optically obscured 
areas. This latter is very dependent on sensor fabrication 
techniques and overall sensor topography design. 

The photon-to-charge carrier conversion efficiency is also 
wavelength dependent. One primary effect for silicon based 
detectors is that for wavelengths greater than 0.7 um, the 
silicon starts to become transparent. Also when the thickness of 
the active region exceeds the pixel dimension, lateral charge 
carrier migration can occur resulting in image spreading (c.f. 
Figure A. 15). 

This spreading will divide the event amount two or more 
pixels resulting in an effective increase in the detection 
threshold. Shorter wavelength and/or special sensor material 

selection (e.g. f Pt Si) can minimize these effects. Relative 
values of the total energy within the 7774 and 8683 angstrom 
lines must be obtained before this spreading effect can be fully 
assessed. Larger pixel dimension is one of the geometric trades 
which can be used to offset the image spreading factor. 

Commerical device quantum efficiencies range from 0.12 to 
0.6 with shorter wavelengths in the region of interest having 
higher values. Some isolated cases have claimed quantum 
efficiencies up to 0.8 for back illuminated devices but this does 
net appear to be cor.sistant with the bulk of the 8683 angstrom 
wavelength devices. 

Typical Goal 

Q- Selected range; 0.35 - 0.45 (0.6) 
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LLfc. 9 1 •/.. SOUI> STATE IMAGE SENSOR POR IMAGE RECORDING 



Horizontal Une-spread function mcafurad at various wavelengths. 
Line of UiumiAaUon with half-wjdth van positioned at the 

center of central pixeL 



Calculated and experimental Unagreed functions at 500, 700, 
and 900 ivn. The calculated points were obtained from a two* dimen- 
sional numerical solution of the dtffanon equation by tbe method of 
eoiyugatc gradients. The bittgrated flux foi cedi pixel is normalised 
lo tbe total flux over the central 9 pixels for both calculated end ex- 
perimental values. The calculated points (open drdcs) are connected 
by a bnc, which semi only as a guide The one ugms error oe dm 
experimental points at low light levels co rres pond* to a relative 
tm y on st of 0.002. Peak illumination level was *50 percent of 
seturation. 


FIGURE A. 15 Image Spreading In Array Detectors. 


Ill 




A. 9. 2 Pixel Dimensions: 
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The selection of pixel size is dependent on a number of 
pararmeters including resolution, aperture, effective F/No. and 
orbit altitude as well as manufacturability. When resolutuion is 
the driver, diffraction limits must be considered. As shown in 
Figure A. 16, the approximate diffraction limited pixel size (L ~ 
?. AR/D) for a F/1 system is around 1.6 um: over a decade below 
the pixel size required for the LITMAP application. Hence LITMAP 
is not diffration limited for global coverage. This figure also 
relates aperture diameter to object resolution for various pixel 
sizes. These factors are related by the pixel angular FOV and is 
given by: 


L » ( R « p ) / ( ( f/No) D ) 


Where all values are in meters and for Figure A. 16 

R » 3.6 10 7 m 

f/NO - f/1 

Each L vs D curve applies to all constant values of 0p/(f/No). 
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Hence an f/1, 40 pm 6 p gives the same curve as an f/2 and 80 um 
pixel. Note that for a specific pixel size, there is a limit on 
the aperture size associated with a given resolution, e.g., 11 km 

imposes an aperture limit of about 5cm for a 15 um pixel. Larger 
pixels permit larger aperture. Figure A. 17 replot3 this latter 
figure giving curves of constant resolution. The SNR 
relationship shows that a larger pixel size (hence larger 
aperture) enhances the system performance. Solar background, 
electronics processing bandwidth and signal dynamic range also 
push for larger pixel size capability to enhance the threshold 
performance of the sensor system. 

In addition to shear bandwidth consideration for a given 
focal plane array, there are also manufacturing considerations. 
Although 800 X 800 arrays have been produced, primarily for the 
astronomy community, they have 15 um pixels and are normally 
.•eadout at sub-MHz rates. Discussions with TI indicates that 



ftp - PIXEL SIZE (pm ) 

FIGURE A. 17 Aperture vs. Pixel Dimensions. 


> 



they would probably not be interested in production of such an 
array with larger pixels, say 20 to 30 um. The physical size 
becomes so large that production yields drop dramatically making 
them prohibatively expensive to produce. Dividing the 800 X 800 
in half, i.e., a 400 X 400 sensor with equal chip frame storage 
makes their production more practical. A number of arrays are 
now in production for commercial TV and scientific applications 
which approach the 300 X 400 and 400 X 500 imaging arrays (plus 
frame storage) with pixel dimensions of 22 to 30 um. 

5 p selected dimension: > 25 um also dependent on other 

system requirements. 

A. 9. 3 Sensor Noise: e 

The primary intrensic sensor noise contr ibut ions are dark 
current and amplifier noise. Analysis indicate that for daytime 
operation, A leakage of 40 NA/cm^ would only contribute less than 
10} more relative to the solar background noise. Existing data 
indicates less than 26 NA/cm^ at 25 C maximum with 8 NA/cm^ being 
typical. Hence, for the short integration times associated with 
LITMAP, there probably will be no cooling requirements due to 
dark current. 

Amplifier noise is a different problem. For sub MHz readout 
rates, as applied in astonomy for example, the preamplifier noise 
is on the order of a few tens of electrons rms* As the readout 
frequency increases, in an attempt to minimize the off focal 
plane electronics required, the preamp noise increases. Special 
on chip preamp designs can be taylored to the desired operating 
frequency regime and special techniques such as corralativn 
double sampling can be applied to minimize the associated readout 
noise. Systems optimization often results in about equal noise 
contributions for the most prominent sources, in this application 
two key contributions are background and readout (prsamp) noise. 
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Desired limits on readout noise will be a contributing factor on 
defining the upper readout frequency for a given array 
configuration . 

The noise factor e was defined as all other noise 
contributions divided by the background contribution (all in 
absolute values, not rms). 

f - Desired value: <1 must be determined for each 

configuration 

A . 9 . ^ In tegration Time: T'i 

The integration time is a delicate trade between (among 
other) excessive noise and power dissipation at higher 
frequencies versus increasing electronics arid required full well 
capacity to hande the solar background at lower bandwidths. 
Figure A. 18 illustrates the relationship between integration 
time, readout rate and pixels per readout port. The partitioning 
given on the right side is based on a total array of 800 X 800 


10 MHZ 


INTERGRATION TIME (MS) 


FIGURE A. 18 Focal Plane Partitioning vs. Sample Time 









(640,000 pixels), A preliminary analysis indicated the desire to 
hold the integration time to a nominal 4 milliseconds. With an 8 
to 10 MHz readout rate, the full array would have to be 
partitioned and provided with 16 output ports. For the four 
array configuration, each array would have four ports. Further 
system optimization may result in 5 MHz readout rate per port, 
hence longer integration time (probably not acceptable) or 
increase partitioning. Electronics weight, power, volume and 
cost must be traded against array manufacturing constraints. 
Four to eight ports per array chip appears to be achievable. 


Analysis discussed later shows that for SNR, r.here is an 
optimum integration time some where around a half to one 
millisecond. From physical considerations, o longer time must be 
utilized . 

A consideration for the upper integration time limit is the 
rate of change of the background. Assuming geosynchronous orbit: 

o Day/Night terminator: approximately 1()3 mi/hr 

or 444 m/sec. Time to move: 


Percent of pixe l FOV 
10 
1 

0. 1 


Seconds Frames 

2.5 625 

0.25 63 

0.025 6 


o Wind (Clouds): less than 10 2 mi/hr 
or 44 M/sec. Time to move: 


Percent of pixel FOV 

Seconds 

Frames 

10 

25 

6250 

1 

2.5 

625 

0.1 

0.25 

63 


Note: vertical cloud growth can occasionally exceed the 

44 M/sec given above. 

For time averaging of the background, a futher restriction 
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may be necessary. For example, for a motion of 0.1% of a pixel 
in ten frames (i.e., 0.1% in MO milliseconds) the terminator 
motion maybe a problem. Satellite jitter must also be assessed 
carefully . 

For LITMAP objectives, the upper integration time will 
probably be more limited by sensor full well capacity rather tha 1 
background changes. 

T^ - Selected baseline : MO X 10"3 seconds 

M - Partitians : 16 

f - Frequency : 10 MHz 

A. 9. 5 Full Well Capacity: FW 

ThxS factor does not explicity appear in the SNR 
relationship but is one of the limiting factors on the upper end 
of the integration time. As background and other noise 
contributions fill the sensor storage area, the allowable signal 
dynamic range decreases. For excessive integration times, the 
ability to handle the signal will go to zero and the pixel 
storage will saturate. 

Some representative full well capacity data is illustrated 
in Figure A. 19 for several commerical ICCD's. The virtual phase 
a"d some three phase devices fall along the lower solid curve 
while the four phase devices fall along the upper curve. 
Although these charge capacities are concept dependent, there is 
some latitude in well depth to adjust the full well capacity. 
Figures A. 20 to A. 22 relates allowable "noise" accumulation to 
signal dynamic range. The signal synamic range is defined as 
maximum signal divided by minimum signal where the mimimum 
detectable signal is defined as: 

N s (peak) s SNR*N n (RMS). 

Figure A.20 illustrates that larger full well capacities permit 
larger allocations to background and noise. The large capacity 
becomes more important as larger dynamic ranges are required 
coupled with larger backgrounds (e.g,, with longer integration 
times). As seen in Figure A.21, this larger capacity allocation 
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2x 10 5 5x 10 5 10 6 


Z AS PERCENT OF FULL WELL IE/FW) 


FIGURE A. 21 Dynamic Range as Percent of Full Well 


to noise is not only a result of an increase in the absolute fu 1 ! 
well capacity, but a larger percent of the full well can be 
allocated to backgroundand noise as the full well capacity 






increases . 
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Finally Figure A. 22 provides curves for fixed dynamic ranges 
and provides the relationship between absolute electron count for 
background (and other noise) as a function of the full well 
capacity. This latter figure illustrates that as the full 
capacity goes up two decades, the allowable background goes up 
2.5 to 4 decades for dynamic ranges of 10 tc 1 0 3 . The 
conclusion: full well capacity is a very important factor when 
dealing with large backgrounds. For a dynamic range of 100, full 
capacities of at least 2 X 105 are almost required with values 
above 5 X 1 0 5 being highly desirable. As shown in Figure A. 19, 
these larger capacities primarily comes from larger pixels with 
some dependency on the CCD technology utilized. The upper curve 
in Figure A. 19 can be represented approximately by 

FW = 1.1 X 103 Ap (um 2 ) 

where Ap is the pixel area. The allowable background can be 
estimated by (c.f. broken curve): 



FIGURE A .22 Absolute Noise Limits Relative to Full Well Capacity. 
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N n 2 = I FW= -7.86 X 10 1 * + 629 A p (um2). 

The lower solid curve in Figure A. 19 is roughly one third of the 
upper curve full well capacity, i.e. for the lower curve 

FW S 0.37 X 10 A p (um 2 ) 

I FW = -2.62 X 1 0 1 * + 210 A p (um 2 ). 

These sets of curves bound most presently existing commercial 
ICCD devices. 

FW - Selected baseline value: > 200,000 electrons 

A. 10 THRESHOLD VERSUS DYNAMIC RANGE 

Whenever the technology is being pushed, th*. inevitable 
question arises as to if performance can be traded versus 
scientific objectives to relieve some technology stressing 
issues. LITMAP is no exception and one of the obvious trades is 
that of threshold setting versus data set completness. 

A previous section discussed observed source strengths. 
This section addresses the relative impact of raising the 
threshold on the loss of data. Figure A. 23 illustrates several 
log-normal distribution curves with the dynamic range being 
defined between the 10% and 90% probability points. The value of 

1 on the ordinate is the nominal setting for threshold to provide 
detection of 90% of the events. 

As the threshold is raised, events are lost at the lower end 
of the curve. A factor of three increase in the threshold for a 
10:1 dynamic range would result in a loss of about 40% more of 
the data. As the dynamic range increases, less sacrifice occurs 
for a factor of two or three increase in the threshold (hence 
SNR). For the LITMAP dynamic range design goal of 100:1, an 
increase of two or three could be tolerated with less than 15% of 
the data being lost. A word of caution is necessary. Present U- 

2 data indicates that the dynamic range may be less than 100 
between the 10% and 90% points, as low as 12:1. That data set is 
very small at present, so care should be exercised when making 
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PROBABILITY 


FIGURE A. 23 Relative Dynamic Range vs. Probability. 


decisions on thresholds. Enough data is available to indicate 
that the above cloud observations are definitely lower than the 
1000:1 dynamic range for below cloud observations. Based on the 
NASA design goals of 100:1 dynamic range the following system 
trade guideline has been assumed where E s is the lowest signal 
that is desired to be detected. 

Ej “ Threshold: Ej < 3hEs 
A. 11 PARAMETER SELECTION SUMMARY 


• A summary of the selected parameter values is given in Table 

j A. 4. These data provided the input for the baseline 
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configuration as given in Table 2.2 of Section 2 and 
foundation of the parametric analysis of Section 3. 


formed the 


TABLE A. 4 Parametric Analysis Selection Summary 


Parameter 

Value 

Units 

Comments 

FOV 

7.37 

Degrees 

Half Angle (50° latitude) 

Pixel 

11 

km 

Sub-satellite (nadir) 

Resolution 

15 

km 

Central CONUS 


20 

km 

50 °N US/CANADA line 

Geometric 

0.6/0.44 


Optimum 

Split, v 

1/0.62 


Subsatellite 


1.35/0.70 


37.5 6 N 


1.82/0.78 


50°N 

\ 

0.54 

Milliseconds 

FWHM 

Pulse 

[T./T. - (T./T,) 2 /3] 

0<T i T l 

Splitting, £ 





[i - (yi^/3] 


T i * \ 


4 

? 

(iJ/m-sr 

Single line. NASA/U2 

X 

8683 

Angstroms 

Alternate 7774 

SNR 

6 


From R(R considerations 

! b 

300 

W/m 2 -pm-sr 

8683A 0 


396 

W/m 2 - M‘"-sr 

7774A® 

Albedo 

0*7 


Both wavelengths 

F/No. 

1.0 



K 

0.4 


Including filter 


5 

Angstroms 



0.15 


Fraction of optical energy 




In single line 

Q 

0.35 


Range 0.35 to 0.6 
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TABLE A. 4 (CONT) 


Parameter 

Value 

Units 

Comments 

6 p 

>25 

microns 

Needs to be large 

c 

< 1 


Device dependent (desire<0.5) 

T i 

4 

Milliseconds 

Integration time 

M 

16 


Total nunfcer of outputs 

F 

10 

** 

Rate per output 

FW 

>2xl0 5 

Electrons 

Full well capacity 

e t 

< 3TJ E s 

nJ/m 2 -sr 

Threshold relative to NASA/U2 data 
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B. PARAMETRIC TRADES 


B . 1 INTRODUCTION 

As mentioned in Section 3, a good feeling for the parameter 
dependencies can be obtained by simple inspection of the basic 
SNR relationship, which in Section 3 was inverted to provide the 
attainable spectral threshold energy for a given set of 
parameters. This inspection approach can be used whenever all 
terms of interest interact in a product/division manner. Hence, 
once a given quantitative value is obtained for the threshold 
energy for a given set of parameters, the adjustment of this 
threshold can be made by inspection. 

This straight forward approach is not as effective when 
addati ve/subtractive factors are involved such as when more than 
one noise term must be considered and for the time and spatial 
terms which have interactive factors. The following sections 
will provide some representative parametric trades which were 
used to determine what parameter set would permit optimization 
of the LITMAP sensor system. The final selection of an 
"optimized" parameter set is constrained by the practical limits 
placed on each parameter as defined in Appendix A. 

B. 2 BASELINE PARAMETER SET 

Appendix A defined a practial parameter set which was deemed 
to be within the limits of the existing technologies and have a 
reasonably low risk for implementation but did not necessarily 
meet all of the requirements of Table 2.1. This parameter set is 
summarized in Table 2.2 as the Baseline. All parametric analysis 
started from thi3 baseline with the intent to determine how the 
overall system (optics/sensor/electronics) could be optimized 
relative to cost, complexity, schedule and risk. 

A graphical presentation of the spectral threshold energy NE 
(J/m 2 -sr) as a function of integration time is shown in Figure 
B. 1 Computer printout formats for the baseline data set is 
illustrated in Figures B.2 (day) and B.3 (night). Note in Figure 
B.1 that the recent NASA/U2 90J threshold value of 4 uJ/m 2 -sr is 
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T NT EGFviT T ON T T ME 1 MOEC • 

8.1 Threshold Energy for 8683A as a Function of Integration 

Time for the Baseline Set of Parameters Given in Table 2.2. 


approached under night time conditions but a value more near 11 
uJ/m^-sr (at 4 milsec) is attainable for worse case daytime 
conditions. From the dynamic range discussion in Appendix A, 
even a factor of 2 or 3 increase in threshold for an event 
dynamic range of 100 would result in loosing less than 10% of the 
data that normally would be attained within the 10 - 90% dynamic 
range band. Previous analysis by others have shown more 
optimistic threshold values due primarily to two factors. First, 
much more optimistic assumptions were taken for parameters such 
as sensor quantum efficiencies than this sturdy had deemed 
acceptable based on commerical devices (c.f. Appendix A for 
discussions). Second, a major contribution to the lower values 
attained in this study was the inclusion of the spatial and time 
splitting factors which split lightning pulses, hence lowering 
their apparent energy. The dashed curves in Figure B.4 result 
when the splitting factors are excluded. These curves are about 
a factor of two below the more realistic curves (solid). A third 
factor which was not included but will have some influence is the 
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LIGHTNING MAF PEP NOISE CALCULATIONS 


F'i ■ el size on Earth 'edge of square 1 

1 1 

i M 

Source size 

1 1 

i hi 

Pi,:el size of CCD 

Tn 

M: crons 

Tel ©scope f /number 

1 


Tel ©scope foe ail length 

9.70 

Crt 

Tel ©scope di ameter 

9. 70 

CM 

Solar lrridiance 

ZOO 

W'M I SF i 

Cloud albedo 

m 7 


Optics transmission, including filter 

. 4 


Wavel ength 

868 . Z 

NM 

Filter bandwidth 

• 5 

NM 

Geometric image splitting factor 

.60 


Lightning Pul sewi dth 

.54 

msec 

Requi red 5/N 

6 


Du an turn ef f 1 c 1 enc y 

.35 


Satellite al ti tude 

35600 

1 M 

Pac) ground 

407300 

el ec t r on s 

Amplifier Noise (rmi) 

150 

Cl ectrcns 



INTEGPhTIOH TIME 'MSEC 


T 

T-FACTOR 

PAO GROUND 

SIGNAL 

PAD' GROUND 

SIGNAL 

TOTAL 



(ELECTRONS) 

(ELECTRONS) 

NOISE 


NOISE 





<RHS EL) 

uJ/M 2SR 

(RMS EL 

.5 

.64 

22687 

1293 

150.62 

7.54 

215.5V 

1 

• 62 

45374 

1561 

213.01 

7. 19 

263.54 

1.5 

.88 

68062 

1823 

260.88 

7.73 

303.95 

2 

.91 

90749 

2037 

301.24 

8.35 

339.53 

2.5 

.92 

113436 

2230 

336.8 

8.96 

3^1.7 


.93 

136124 

2407 

368.95 

9.56 

401 . 28 

3.5 

.94 

156811 

2572 

398.51 

10. 12 

428.81 

4 

.95 

181498 

2728 

426. 02 

10.66 

454.67 

4.5 

.95 

204186 

2B74 

451.86 

11.17 

479. 12 

5 

.96 

226873 

3014 

476.31 

11.67 

502.38 

5.5 

.96 

249561 

3147 

499.56 

12. 14 

524.6 

6 

.96 

272248 

3275 

521.77 

12.6 

545.91 

6.5 

.97 

294935 

3398 

543.07 

13.04 

566.42 

7 

.97 

317623 

3517 

563.58 

13.47 

586.2 


FIGURE B.2 Computer Printout Format for Basel Ine/Day 
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Ri*el n:e on Earth (edge of square) 

1 1 

t m 

Sourct lire 

1 1 

I m 

Ri::el siz« of CCD 

30 

Ml cron?. 

Telescope f /number 

1 


Telescope ^ oca 1 1 eng t h 

9. 7v 

CM 

Telescope diameter 

9. 7 <; 

CM 

Solar irradiance 

0 

W/M I SR i 

Cloud albedo 

.7 


Optica t-anami ssi on , including -filter 

.4 


Wavelength 

86e. 7 

NM 

Filter bandwidth 

. * 

NM 

Geometric image splitting factor 

. 62 


Lightning Pulsewidth 

. 54 

msec 

Required 5/N 

6 


Quantum efficiency 



Satellite altitude 

35600 

y n 

Bad ground 

487500 

el fctroni 

Amplifier Noise <rms) 

150 

Electrons 



INTEGRATION TIME -MSEC- 


T 

T-FACTOR 

BAOr GROUND 

SIGNAL 

BACK GROUND 

SIGNAL 

total 



(ELECTRONS) 

(ELECTRONS) 

NOISE 


NOISE 





(RMS EL) 

uJ/M 2SR 

(RMS EL) 

.5 

• 64 

0 

91 B 

0 

5.35 

153.02 

1 

• 82 

0 

918 

0 

4. 17 

153.02 

1.5 

• 88 

0 

918 

0 

3.89 

133.02 

2 

.91 

0 

918 

0 

3.76 

133.02 

2.5 

.92 

0 

918 

0 

3.69 

133.02 

•v 

.93 

0 

918 

0 

3.64 

133.02 

3.5 

.94 

0 

918 

0 

3.61 

133.02 

4 

.95 

0 

9J8 

0 

3.58 

133.02 

4.5 

.95 

0 

*18 

0 

3.56 

133.02 

5 

.96 

0 

918 

0 

3.55 

133.02 

5.5 

.96 

0 

918 

0 

3.54 

153.02 

6 

.96 

0 

918 

0 

3.33 

133. 02 

6.5 

.97 

0 

918 

0 

3.32 

133.02 

7 

.97 

0 

918 

0 

3.31 

133.02 


FIGURE B.3 Computer Printout for Baseline/Night 
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FIGURE 8.4 Baseline with ( — ) and without (— ) Spatial 
and Time Splitting Factors 


frame transfer dead time which will probably be between 10X and 
3 0 X depending on the selected integration time. Due to 
electronic/sensor system constraints, the 10% level is being 
sought. These r actors just emphasis that parameter evaluation 
and selection must be done carefully. The values selected permit 
a certain amount of latitude for improvement in some areas to 
compensate for some inevitable degradation of other parameters. 

o Lightning event splitting (time/spatiai) is important. 
B. 3 OPTIMISTIC SET 

Although many parameters can be adjusted slightly, three 
parameters have been identified which, although will result in 
some increase in risk and possibly cost, would provide some 
enhancement of the LITMAP performance. These parameter (K, Q, and 
pixel size) adjustments are provided in Table 2.2. Figure 8.5 
shows that these changes will reduce the attainable threshold 
value to near 5 uJ/m2-sr. Hence desired performance based on the 
given assumptions, including source image size equal to pixel 
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FIG 1 IRE B.5 Optimistic Performance 


size, approaches the desired detection threshold capability. 

o Increased performance can be attained with some 
increase in risk and possibly cost. 

B. 4 LIGHTNING PULSE DURATION 

The effect of dividing lightning events between two frames 
is illustrated in Figure B.6 along with effects of pulse 
duration. All of the curves in this appendix have this general 
shape. A finer step and larger range for the integration time 
was used in Figure B.6 to emphasis the shape of this curve. From 
hardware contraints, times shorter than one millisecond would not 
be practical so mest of the parameter trades were done with fewer 
data points and longer times. 

The mimimum threshold point is at the minimum of the curves 
in Figure B.6. The left hand (small integration time) rise is 
due to the time splitting factor and the gradual right hand rise 
is due to the longer integretion of the solar background. The 
minimum occurs where the integration time equals the pulse 
duration time. To avoid large changes in threshold as a 
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FIGURE B.6 Pulse Duration Effects 


function of pulse duration, integration times longer (rather than 
smaller) than the pulse duration is desired. This integration 
time requirement is also driven by hardware, hence the analysis 
in Appendix A resulted in a baseline selection of 4 milliseconds 
as the integration time. This selection somewhat minimizes the 
variations due to pulse durations from sub milliseconds to one or 
two milliseconds. Not that the longer pulses will have some 
increase in their threshold, but these pulses tend to be the more 
intense pulses, hence will be detected inspite of this increase. 

o Pulse duration and time splitting considerations 
drive the integration time toward several milliseconds. 
This result is also compatible with some of the 
hardware constraints (e.g. data rate). 

B. 5 FILL FACTOR 

The source size relative to individual pixel FOV dimension 
influences the threshold performance capability of the LITMAP 
instrument. Figure 3.7 illustrates the effects where all 
parameters were set at their baseline values (including fixed 
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FIGURE B.7 FOV Effects - Day 


source size) and the FOV varied. The 0°N (subsatellite point) 
and 50°N (maximum viewing angle) limits are included. The 
analysis was done such the maximum permissible aperture was 
always used as constrained by the FOV and detector pixel size. 
This gave the minimum threshold for a given set of conditions. 
As expected, the FOV to source size has strong effect on the 
threshold capability. It is thus emparative that the effective 
source size be quantified. This will be discussed a little more 
in a subsequent section. 

For the present situation, with a fixed effective source 
size, the background area viewed is larger as the viewing angle 
moves away from the sub-satellite point. A threshold degradation 
of about ?.7 results at the maximum viewing angle. 

Spatial splitting has been included and does have an effect 
which is very apparent in the nighttime plots of Figure B.8 
corresponding to the same conditions of Figure B.7. Without the 
splitting factor the nighttime capability would be only amplifier 
noise dependent, hence constant (without noon light etc.). But 
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FIGURE B.8 FOV Effects - Night 


as seen in Figure B.8, the nighttime capability degrades 
(threshold increases) as the FOV increases, primarily because of 
the effect through the spatial splitting factor which requires 
larger effective source strengths as the FOV increases relative 
to source size. 

The composite result for the FOV effect coupled with the 
spatial splitting falls between a linear and square dependancy as 
illustrated in Figure B.9 where the solid line is for computed 
values. All curves were made equal at the nominal value of 11 km 
baseline sub-satellite resolution point. 

o Off nadir viewing angles will have a threshold 
degradation of up to a factor of about 3. 

B.6 PIXEL SIZE 

One of the key parameters that can be varied to enhance the 
L1TMAP instrument performance is that of sensor pixel size. 
Larger pixels permits larger apertures and full well capacities. 
This latitude must then be traded with achievable array size 
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which usually results in building a mosaic of smaller arrays. 
This mosaic approach introduces its own complications including 
optical/focal plane complexities. 

Performance for pixel sizes from 15 to 40 microns are 
illustrated in Figure B. 1 0 for daytime operation. The 15 micron 
pixel size as used for the large 800 X 800 NASA/JPL Galileo 
project is not suitable for LITMAP, both from sensitivity 
(limited aperture size) and full well (background handling) 
requirements. Figure B. 1 1 again illustrates that the larger 
aperture permitted by a larger pixel results in a lower 
detectable energy level. 

o Large pixel size is necessary (> 20 microns) 

B. 7 READOUT NOISE 

The process of extracting the photon produced electrons 
from the sensor array will induced noise contributions in 
addition to those associated with solar background and signal. 
The primary contributor, as indicated in Appendix A, is the 
readout amplifier associated noise. This noise contribution is 
often specified as the number of rms noise equivalent electrons. 

For the baseline conditions, the effect of varying the 
amplifier noise is smmarized in Figures B.12 (day) and B. 1 3 
(night). For normal solar background and cloud albedo, a 
significant contribution from the amplifier noise does not occur 
until somewhere above 200 e rms. From Figure B.12, the present 
baseline concept could accommodate up two or three hundred 
electron rms noise contribution without serious degradation of 
performance. But as noted in Figure B. 1 3 » as the amplifier noise 
increases, the night performance degrades. As this contribution 
approaches 400 e rms, the day/night performances will be about 
equal for a 4 millisecond integration time. For larger noise 
contributions, the performance will become amplifier limited, not 
solar background limited. 

o Amplifier noise must be kept significantly below the 
solar contribution if enhanced night performance is 
desired (i.e., N amp (rms) < N sol (rms)/2). 
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FIGURE B.10 Pixel Size Effects - Day 



FIGURE B.ll Pixel Size Effects - Night 
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FIGURE B.12 Amplifier Noise Effects - Day 
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FIGURE B.13 Amplifier Noise Effects - Night 


* 

t 

# 

4 


137 




B.8 WAVELENGTH 
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Primary focus has been on the 8683A wavelength. Other 
candidates exist and should possibly be evaluated. The baseline 
performance at 8683 and alternate 7774 angstroms were compared. 
The following values were assumed: 


Parameter 

Baseline 

Alternate 

Units 

Wavelength 

8683 

7774 

A 

Q 

0.35 

0.5 

- 

IB 

300 

396 

w/m2-sr-um 


Figure B.14 illustrates that the threshold values are about the 
same for the conditions considered. Other factors such as sensor 
materials and/or larger absolute strength for one of the lines 
may result in a preference. At the present, the 8683 A line has 
been baselined although it has a somewhat lower detection quantum 
efficiency when silicon dioxide is used for the structural 
elements. If tin oxide is used, as for the Westinghouse devices, 
both spatial line detection efficiencies are about equal and 
fairly high (> 50$). Hence the tin oxide material would be a 
good candiate for the "optimistic” system. 

o Relative optical magnitudes of the spectral lines are 
needed . 
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FIGURE 8.14 Wavelength Selection 
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B.9 FILTER BANDPASS 
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The filter bandpass must be as narrow as possible to reject 
the solar background. Birefringent filters have been developed 
that can provide large FOV coverage while maintaining narrow 
passbands but suffer from lower transmission, high cost and high 
risk. Interference filters, although more technologically 
available, exhibit a FOV angular dependency which does not 
accommodate a large FOV f18°) and narrow bandwidth (<5 A) 
simultaneously. 

The effect of optical bandpass or. the daytime threshold 
performance is illustrated in Figure B.15. A 25A bandpass would 
degrade the system performance by a factor of 2 relative to a 5A 
filter. This degredation is in addition to the factor of 2 which 
occurs if optical splitting factors are not included. 

As a result of these considerations, the FOV was divided 
into smaller portions which permitted the use of interference 
filters possessing around 5A passband to be used. If the' need 
is sufficient, some improvement could be realized through the use 
of a birefringent type of filter. 

o Narrow bandpass is necessary (<5A). 
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FIGURE B.15 Optical Bandpass Effects 
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B. 10 EVENT SOURCE SIZE 
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Physical, complexity and cost constraints place restrictions 
on the allowable number of array elements that can be utilized. 
This limited array size coupled with full disc coverage places a 
severe restriction on the source energy that can be resolved 
during day operations. As discussed earlier, this detection is 
also complicated by the geometric splitting. 


The baseline concept assumes that the source and pixel FOV 
are of equal size, 11km. Figure B.16 illustrates the increase in 
the required threshold energy as the source moves from 11 km 
dimensions (lower curve) to 5 km. The night performance is not 
as bad as illustrated in Figure B. 1 7. If the 4 uJ/m^-sr is a 
true indication of the source strength, and if the effective 
source size is considerably less than 11 km square, then the 
threshold will be much larger. This condition would result in 
pushing the LITMAP instrument design hard toward the technology 
limits tc provide the performance desired from a geosynchronous 
orbit. An alternative would be to decrease the pixel FOV to 
match source size. If the full disc coverage then requires too 
many array elements, a move to a very LITMAP compatible 
step/stare concept could be made. 


o Event source size must be quantified. 
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FIGURE B . 16 Source Size Effects - Day 
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FIGURE B.17 Source Size Effects - Night 


B. 1 1 EXISTING DEVICES 

The capability of existing imaging array devices were 
assessed to put the overall system trade evaluations into 
perspective with the real world. A 22.4 micron TI (TV: IC-202) 
device is compared with a RCA (SID 501D) 30 micron device. As 
can be seen in Figure B . 1 8 , the RCA device comes close to the 
LITMAP baseline limit. The smaller pixel device has a higher 
threshold capability. Performance of other devices are 
illustrated in Figure B. 1 9 (8683A) and Figure B.20 (7774A), For 
ease of wavelength dependency comparison, Figure 21, 22 and 23 
illustrates the 8683A and 7774A results for RCA, Westinghouse and 
Hughes. Where the detection quantum efficiencies are nearly 
equal for the two wavelengths (as in the Westinghouse device) the 
8683A line comes out ahead, otherwise the 7774A line has a small 
edge. The device pixel sizes and quantum efficiencies are given 
in a device characteristics table located in Section 5. 

o Existing commercial devices approach the LITMAP 
baseline performance specifications (array 
partitioning required to handle data rate). 
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FIGURE B- 18 Actual Devices - RCA, TI 




FIGURE B.19 Actual Devices - 8683A 
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FIGURE B.?0 Actual Devices - 7774A 
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TABLE B.21 RCA: 8683A and 7774A (SID 5010) 
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FIGURE B.22 Westinghouse: 8683A and 7774A (5040) 
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FI6URC B.23 Hughts: 8683A and 7774A (H4068) 
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B. 1 2 SUMMARY 
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This section has illustrated the threshold dependency on the 
various sensor design parameters. These analyses were also 
weighed against electronic and physical constraints as well as 
cost and schedule. Table B,1 summarized the results of this 
append ix . 


TABLE B.l Parametric Trade Results 


FACTOR 


COMENT 


o Baseline 
o Optimistic 

o Pulse Duration 

o Fill Factor 

o Pixel Size 

o Readout Noise 

o Wavelength 


Lightning event splitting (time/spatial) is important. 

Increased performance can be attained with some increase 
in risk and possible cost. 

Pulse duration and time splitting consicsrations drive 
the integration time toward several milliseconds. This 
result is also compatible with some of the hardware 
constraints (e.g, data rate). 

Off nadir viewing angles will have a threshold degradation 
of up to a factor of about 3. 

Large Pixel size Is necessary (> 20 microns) 

Amplifier noise must be kept significantly below the 
solar contribution if enhanced night performance Is 
desired (i.e., N a|np$ (rms) <N $ol (rm$)/2). 

Relative optical magnitudes of the spectral lines are 
needed. 


o Filter Bandpass Narrow bandpass is necessary (<5A). 

o Event Source Size Event source size must be quantified. 

o Existing Device Existing commercial devices approach the LITMAP 
baseline performance specifications (arrav 
partitioning required to handle data rate). 


S. GOVERNMENT PRINTING OFFICE 1983-746-070/81 


145 


